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ABSTRACT 

The dynamics associated with the perception of orienta- 
tion. were modelled for near-threshold and suprathreshold 
vestibular stimuli • A model of the information available 
at the peripheral sensors which was consistent v?ith available 
neurophysiologic data was developed and served as the basis 
for the models of the perceptual responses . As a preliminary 
assumption the central processor was to utilise the 

information from the peripheral oonsors in an optimal (mini- 
mum moan ntjuaro. error) manner to produea th@ perooptuai esti- 
mates of dynamic orientation. Thi© assumption, coupled with 
the models of sensory information, d®t©rmin@d the form of 
the model for the central processor. Comparison of model 
responses with data from psychophysical experiments indicated 
that while little or no central processing may be occuring 
for simple suprathreshold canal stimulation, a significant 
portion of the dynamic response to translational accelerations 
must be attributed to the central processing of otolith infor- 
mation. 

The fundamental mechanism which underlies the phenomenon 
of vestibular thresholds was studied experimentally by 
testing the response of subjects to a near threshold stimulus 
consisting of a velocity step-ramp proportional to the sum of 
the subject’s velocity step and acceleration step thresholds. 
Experimental results indicated that canal thresholds could be 
accounted for by a model of central processing consisting only 
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of: an optimal processing of afferent fifing rates in additive 
noise with no necessity for peripheral dead zone nonlinearities. 
Quantitative models of threshold detection were developed which 
correctly predicted threshold levels (75% correct detection} and 
response latencies for rotational stimuli. It was found that 
the same detector could be used to model the threshold responses 
resulting from translational stimuli. 

The illusions of static orientation were studied and it 
was shown that they were consistent with a simple vector 
transformation which could be associated with differences in 
the processing of signals arising from stimuli in and stimuli 
perpendicular to the "utricle plane". A model was developed 
which incorporated this difference and which was capable of 
predicting the perception of orientation in an arbitrary 
static specific force environment. 

The problem of integrating information from the semi- 
circular canals and the otoliths to predict the perceptual 
response to motions which stimulate both organs was studied. 

A model was developed which was shown to be useful in pre- 
dicting the perceptual response to multi-sensory stimuli. 
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Chapter I 
INTRODUCTION 

The research effort which is described in this thesis 
was undertaken to increase our understanding of the • phenom- 
enology associated with vestibular perception. More 
specifically , this thesis attempts to separate the processes 
which underlie this perception into two cascaded but funda- 
mentally distinct elements, namely: 

1. the peripheral sensors and the associated neural 
processes which determine the afferent response 
to external stimuli; and 

2. the processing by the higher centers of the 
information available from the first order afferent 
responses . 

These elements' can be considered as both physically separable 
and, in terms of our methodological approach to modelling 
them, philosophically separable. In the case of the first 
element the mechanical dynamics of the sensors can be and 
are distinguished from the dynamic effects associated with 
the first order afferent processes. In the case of the 
second element a distinction is made between the threshold 
stimuli which must be processed by a detector and supra- 
threshqld stimuli which must be processed by an estimator. 
Supra threshold stimuli are further divided into those which 
involve the integration of more than one sensory modality 
and those which do not. The motivation for this research 



and the methodology used to approach the fundamental problems 
Involved in the modelling of vestibular perception are 

discussed in this chapter followed by a brief introduction 
to the contents of the remaining chapters . 

1 o 1 Motivation for Research 

The products of modern technology, especially those 
associated with the advancements made in aerospace vehicles, 
have engendered a rapid increase in the need to understand 
man's reactions to motion environments which are completely 
alien to his pre-twentieth century experiences. Obvious 
examples of such environments are the prolonged zero-g 
environments made possible by space vehicles and the rapidly 
varying high g environments of modern military aircraft and 
rocket launch vehicles* Less obvious but equally important 
are man's reactions to the motions arising from commercial 
and general aviation aircraft, ships, tall buildings (sway) 
and motion based simulators. Since the vestibular system is 
man's primary non-visual inertial orientation sensor, its 
central importance to any understanding of man’s capability 
to function effectively in these motion environments is 
clear. 

If a model of stimulus detection is developed 

for vestibular perception which is capable of giving reasonable 
estimates of the detection probabilities as a function of 
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time; for arbitrary near threshold stimuli then predictions 
can be made which have significant importance for several 
seemingly unrelated problems . Among the most prominent of 
these are the following: 

1. How does one maximize the fidelity of a motion- 
based simulator while minimizing the requirements 
for translational motion so that simulator costs can 
be reduced? While the techniques involved in 
optimally utilizing a given amount of lateral 
motion capability can become quite sophisticated 

a thorough understanding of the dynamics of threshold 
perception is necessary if maximum fidelity is to be 
achieved. 

and ■' 

2. What are the constraints which must be placed on 
the structural motions of tall buildings to insure 
the comfort of the buildings' occupants? 'fhis question 
is of groat importance in the design of tall buildings 
since that design (and therefore the construction 
costs) will bo very sensitive to the constraints 
imposed, Since these motions are typically quite 
small it is necessary to have a general model for 

the detection of near threshold motions if reasonable 
trade offs are to be made. 

The need for a model to predict the subjective perception 
of dynamic orientation for suprathreshold stimuli which 
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involve stimulation of both the semicircular canals and the 
otoliths is even greater- Such a model could be used to 
study postural control , to evaluate the ride quality of a 
wide variety of transportation vehicles , to predict the 
reactions of pilots during unusual maneuvers/ to predict 
the incidence of motion sickness, and to evaluate many of 
the illusions of motion or orientation which arise due to 
unusual g forces or sustained rotations. In addition, the 
development of such a model draws upon and may contribute 
to knowledge of sensory/neural physiology . For example 
investigators who are studying the neural processing centers 
of the brain may find interesting parallels between the 
interactions of sensory information they discover and the 
mathematical transformations required by the model. If 
the time comes that a one-to-one correspondence can be made 
between a mathematical model of man's perceptual responses 
and the processes seen in the brain then it might be possible 
to predict the site of neurological disorders based upon 
the response of patients to controlled stimuli. 

1.2. Approach to the Problem of Vestibular Modelling 

Since this research is concerned with the processing 
of information which is relevant to the perception of 
dynamic orientation it is important to elucidate carefully 
exactly what information is being considered. While the 
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problem of integrating the signals available from every 
sensory system which provides information pertinent to the 
perception of dynamic orientation is an important one, 
it is a task which, due to its extremely broad scope must 
await the solution of more restricted problems. The models 
developed in this thesis exclude any information gained 
from non-vestibular sensors during the time of stimulus 
exposure. Specifically excluded are visual, tactile, pro- 
prioceptive, kinesthetic and aural information. Information 
gained prior to the stimulus exposure is considered a priori 
information and in most cases can be handled by the general 
mathematical framework of the models if care is taken to 
account fully for' the nature of that information and how it 
is affected by any pre-stimulus instructions. In addition 
to ii priori information the higher centers have at their 
disposal two other types of information. The first of these 
is the afferent signal available from the vestibular sensors 
which must be processed to obtain the perceptual estimates of 
orientation. To process these sensory signals and to mix 
them optimally with the a priori information requires some 
knowledge of the processes which give rise to the afferent 
signals. It is this knowledge (which includes an internal 
model of both the sensory dynamics and the measurement noise 
processes) that completes the information base available to 
the higher centers. 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 



.( n the proccMin of developing models for the information 
available from the vestibular sensors and for the processing 
of that information by the higher centers a number of issues 
arise which can be dealt with as methodological problems . 

Several of these issues will be dealt with here as illustra- 
tions of the approach to perceptual modelling taken in this 
research . 

One question which arises immediately in any modelling 
effort concerns the criteria which will be used to select 
the form of the model and its parameters. The answer to this 
question depends on the amount of .knowledge available about 
the. physical system being modelled. Since a significant 
amount of qualitative and quantitative knowledge is available 
concerning the mechanical and afferent dynamics of the 
vestibular sensors this information will be used as much as 
possible in modelling their dynamic response. On the other 
hand the knowledge available about the internal structure 
and organization of the central processor is much more 
limited, so much so in fact that any attempt to develop a 
viable model of subjective perception as a function of afferent 
responses based upon the known neurophysiological structure 
of the brain would probably be fruitless. Upon what information 
then can a model of the central processor be based? In any 
modelling effort in which the physical structure is considered 
completely unknown, a "black box" approach is taken to produce 
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a model which is consistent with the known responses of the 
system to selected input3 . In this case it is reasonable 
to presume that the neural net which forms the central 
processor has evolved to use the available sensory informa- 
tion in a roughly optimal way. This assumption of optimality 
serves to suggest the form of the processor and eliminates 
any problem of non uniqueness (the fact that more than one 
processor might have been capable of producing the required 
predictions) . The model which proceeds from this assumption 
of optimality must then be checked against the known output 
of the system — namely the subjective responses determined 
from psychophysical experiments. 

Another problem which arises involves the issue of 
efferent signals from the higher centers which may dramatically 
alter the afferent response. The mechanism by which efferent 
discharges affect afferent responses is, as of yet, unknown. 
Whatever that mechanism, if the transformation of the 
afferent signal has a unique inverse then the point of view 
taken in this research is that the total information available 
to the higher centers concerning the dynamic orientation 
of the head is unchanged since the' higher centers initiated 
the efferent signal and are aware of its effect on the 
afferent response. Since the goal of this research is to 
develop a model of the available afferent vestibular informa- 
tion and its subsequent optimal processing by the higher 
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centers any invertible alteration of the afferent information 

s 

will not affect the models predictions and therefore can be 
ignored in the development of the model. 

Finally the problem of active versus passive head move- 
ments is an issue which deserves attention. Since an active 
movement of the head is initiated by the higher centers, 
the resultant motion (to the extent that it can be predicted 
open loop by the higher centers)' is available as a priori 
knowledge. Therefore the equivalent information available 
from the peripheral sensors is redundant and adds nothing 
to the a priori knowledge available to the higher centers. 
Since the open loop estimate of head motion is bound to 
contain some errors the afferent response from the peripheral 
sensors should be used to check these a priori estimates. One 
way to accomplish this is to take the difference between 
the expected sensory response (based upon the ' open loop 
estimate of motion and the internal model of the sensory 
dynamics) and the actual sensory response. This is one inter- 
pretation of the corollary efferent discharge or “efferent copy. 
The resulting signal can then be processed to estimate the 
errors associated with the movement. The models developed 
in this thesis should therefore only be used in an error 
correcting mode for motions which are initiated by the 
higher centers. 
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Additional problems similar to those described above 
arc dealt with in a similar fashion when they arise in the 
course of developing the models . 

1.3. Thesis Organization 

Chapter Two summarizes the structure, function and 
orientation of the semicircular canals and the otoliths. 

Chapter Three derives a model of the information available 
at the first order afferent level of both the semicircular 
canal system and the otoliths. These models are then coupled 
with optimal estimators to yield predictions of subjective 
perception for simple noninteracting stimuli. 

Chapters Four, Five, and Six develop models for the 
detection of near threshold stimuli. Chapter Four describes 
an experiment which was conducted to determine the fundamental 
mechanism underlying the threshold phenomenon. Chapters Five 
and .Six develop quantitative models for the detection proces- 
ses associated with rotational and translational motions 
respectively. 

The seventh chapter investigates the illusions of static 
orientation as a function of body position and the strength 
of the gravito-inertial field. A simple mechanism is pro- 
posed which accurately predicts these illusions and indicates 
that they most likely have a common origin. 

The problem of integrating the information from both 
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the semicircular canals and the otoliths to arrive at a 
single set of preceptual responses is investigated in Chapter 
Eight. A model of sensory integration, is proposed which can 
be used either qualitatively or quantitatively to predict 
the subjective perceptions associated with interacting 
stimuli. 

Finally^ Chapter Nine summarizes the conclusions which 
can be drawn from this research and suggests possibilities 
for further experimental and analytical investigation. 
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Chapter II 

THE HUMAN VESTIBULAR SYSTEM 

The purpose of this chapter is to introduce the reader 
who is unfaniliar with the vestibular sensors to the 
basic structural organization and physiologic function 
of these organs. A more in-depth introduction is available 
in references 9,61 and 80. 

2.1. Semicircular Canal System 

The semicircular canals are the primary nonvisual sensors 
of rotational motion with respect to inertial space. They 
consist of three approximately circular toroidal canals 
whose axes form a roughly orthogonal set. The membranous 
canals are suspended in a fluid (perilymph) in the temporal 
bone of the skull adjacent to the auditory portion of the 
inner car. Figure 2.1 illustrates the entire inner ear (in- 
cluding vestibular and auditory portions) and Figure 2.2 
indicates the orientation of the canals relative to the head. 
The semicircular canals are filled with a water- like fluid 
called endolymph which, due to its inertia, tends to lag 
behind the motion of the canal walls when the head undergoes 
angular acceleration. When the endolymph moves relative to the 
canal it tends to displace the cupula which obstructs an 
expanded section of the canal called the ampulla. This dis- 
placement of the cupula is detected by sensory hair cells at 
tnc ease of tne cupula which in turn produce a change in 
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Fi yurt/ 2.1 Diagram of Hainan Inner liar 

(Abbott Laoorutories R<af . 1) 
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KEY 

H RIGHT 
L LEFT 

H HORIZONTAL 
P POSTERIOR 
S SUFERIOB 
(X+) POSITIVE AC- 
CELERATION ABOUT 
THE X AXIS INCREASES 
THE AFFERENT FIRING 
RATE 

(Y-) POSITIVE AC- 
CELERATION ABOUT THE Y 
AXIS DECREASES THE AF- 
FERENT FIRING RATE 


Figure 2,2 Orientations nnrt Polarizations of 
the Semicircular Canals 
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the firing frequency of the first order afferents which 
provide information to the central nervous system- Figure 2-3 
illustrates this process for the horizontal canal- All of • 
the hair cells associated with a particular canal have the 
same polarization, i.e., displacement of the cupula due 
to endolymph flow in one direction will either excite all 
of the sensory hair cells or inhibit them all. 

Since the canals on the right side are essentially 
coplanar with the canals on the left side they are pairwise 
sensitive to angular accelerations about the same axes. 
Investigation of the afferent responses of these sensors 
indicates that a pair of canals which are sensitive to acceler- 
ation about the same axis (e.g. the right posterior canal and 
the loft superior canal) have opposite sensitivities (see 
Figure 2.2B) , so it is presumed that the higher centers 
respond to the difference of their responses. The detailed 
dynamic response of the afferent firing of the semicircular 
canals to an angular acceleration of the head is discussed 
in section 3.1. 

2.2. Otoli th System 

In addition to the semicircular canals, the nonauditory 
portion of each inner ear contains two otolith organs which 
are sensitive to changes in the gravitoinertial reaction force 
(referred to here as specific force). .The approximate 

EEPRODUCIBILnY OF THE 
tiRIOSNAfe MGR IS POOR 
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locution of each of these organs, known as the utricular 
otolith and the saccular otolith, is shown in Figure 2.1. 

Each otolith organ contains a gelatinous layer interspersed 
with calcium carbonate crystals and supported by a large 
number of sensory hair cells. Since the calcium carbonate 
crystals (known as otoconia) have a higher specific gravity 
than the surrounding fluid (endolymph) an appropriate 
acceleration of the head will tend to shift the 
otoconia relative to the bed of sensory cells (known as 
the macula) . When this shifting motion occurs the sensory 
hairs are bent and the afferent fibers which innerviate these 
hair cells change their firing rate. Figure 2„4 illustrates 
the basic structure of the otoliths. 

Motion of the otoconia parallel to the bed of. sensory 
hairs (in Figure 2.4: motion right and left or into and out 

of the page) is normally assumed to be the effective agent 
in eliciting a change in afferent firing. The utricles are 
oriented such that the major plane of their sensitivity is 
parallel to the plane of the. horizontal semicircular canals. 
The saccular organs are oriented so that their plane of 
sensitivity is perpendicular to the horizontal canals (and 
therefore the utricles) and roughly parallel to the 
median plane. Figure 2.5 illustrates the approximate orienta- 
tion of the otolith organs. Unlike the semicircular canals, 
the hair cells in the otolith organs do not all have the 
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AFFERENT NERVE 
FIBERS 


Figure 2.4 Cross Section of Otolith Organ 
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same directional polarization. Figure 2.6 illustrates 
the general morphological distribution of hair cell polari- 
zations for the utricle and saccule. The directional 
distribution of polarizations in the utricle is reasonably 
uniform and thus the utricular otolith can be considered 
approximately equally sensitive to shear forces in any 
direction in the utricular plane. The distribution of 
polarizations in the saccule is much more restricted , with 
the major axis of sensitivity roughly perpendicular to 
the utricular plane. Therefore the saccule can be considered 
as an accelerometer sensitive to changes in the specific force 
perpendicular to the average plane of the utricles. 

The detailed dynamic response of the otolith afferents 
is discussed in section 3.2. 
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POSTERIOR 




Figure 2,6 Morphological Polarization Maps for the 

Saccule and Utricle. of the Squirrel Monkey 
(After Lindemnn, Ref, 42; 
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Chapter III 

MODELLING OF FIRST- ORDER AFFERENT S AND RESPONSE TO 
NON INTERACTING - SUPRATHRBSHOLD STIMULI 

The purpose of this chapter is to develop models of the 
sensory information available at the first order afferent 
level for both the semicircular canal and the otolith systems. 
The resulting models will serve as the informational link be- 
tween the true motion with respect to inertial space and 
higher processing centers in the brain. Since we take the 
view that the higher centers have most likely evolved as op- 
timal or near optimal processors of this information, a 
specification of the relevant sensory dynamics plays a major 
role in determining the overall dynamics of the subjective 
response to vestibular stimulation. Once the first order 
afferent response is modelled for each of the vestibular sen- 
sors and these models are coupled with reasonable models of 
process noise (which represent the a priori information con- 
cerning the statistical nature of the expected input) and 
measurement noise, then the optimal processor can be formu- 
lated and predictions made concerning the subjective response 
to noninteracting suprathreshold stimuli. The phrase “non- 
interacting suprathreshoid stimuli" refers to any supra- 
threshold stimuli which involves no change in the orientation 
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of the subject with respect to the gravitational vertical 
and for which the subject is consciously aware that no such 
change will take place. Generally, this means rotations 
which are performed about an axis parallel to the local g 
vector and accelerations which are performed in a device 
which the subject knows is incapable of rotations out of 
the vertical. The reasons for these limitations will be- 
come clear when stimuli not meeting this description are 
considered in Chapter Eight, 
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3.1 Semicircular canals 

The specification of afferent dynamics for the semicir- 
cular canal system is divisible into several parts. The 
first involves a modelling of the mechanical movement of the 
cupula within the ampullary lumen. The second part of the 
specification concerns itself with how this mechanical move- 
ment is reflected in the neural firing rate in the first 
afferent nerve. Finally, an assessment must be made of that 
portion of the afferent signal which is found to be indepen- 
dent of the stimulus input and which therefore is considered 
to be measurement noise in the context of this modelling 
effort. 

3.1.1 Dynamic Response of Cupula 

The structure and fundamental mechanical operation of 
the semicircular canal is described in section 2.1 and illus- 
trated in Figure 2.3. Two forces act to accelerate the endo- 
lymphatic fluid (which fills the canals) with respect to 
inertial space. The first of these is a viscous drag which 
is proportional to the rate of movement of the fluid with 
respect to the walls of the canal. The second is presumed 
to be a linear elastic restoring force which arises from the 
spring-like tendency of the cupula and/or the membranous 
canal to maintain its resting position or shape, 

REPRODUCIBILITY OF THE 
ORIGINAL, PAGE IS POOR 
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If it is assumed that for normal physiologic motions 
no endolymphatic fluid is allowed to leak between the cupula 
and the inner wall of the ampulla, then the amount of move- 
ment of the cupula is proportional to the angular motion of 
the fluid with respect to the canal walls. This assumption 
seems to be warranted on the basis of observations by 
Steinhausen (Ref. 66 ) and injection micrographs performed 
by Groen, Lowenstein and Vendrik (Ref . 32 ) . Using this 

assumption, the motion of the endolymph relative to the canal 
can be expressed as follows: 


M (0 + 0 .) ~ -V0 - K0 

ec ci' ec ec 


(3.1) 


where 0 


ec 


0 


ci 


M 

V 

K 


angular deflection of the endolymph with re- 
spect to the canal 

angular position of the canal with respect to 
inertial space about an axis normal to the 
plane of the canal 
moment of inertia of the endolymph 
coefficient of viscous drag- 

coefficient of linear restoring force due to 
displacement of the fluid within the canal. 


Equation 3.1 is referred to as the torsion pendulum 
model and was first developed (using different conventions) 
by Steinhausen (Ref. 67 ) after observing the motion of the 
cupula in the pike. The equation is arranged to illustrate 
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that the acceleration of the endo lymph with respect to the 
inertial space (0 _ + 0 . ) is due to the sun of a viscous 

wC Cl 

drag force (~VG) ) and an elastic restoring force (-K0 ) . 

ec sc 

Equation 3,1 is time invariant and can be Laplace transformed 
to yield the following transfer function relating endolymph 
displacement to acceleration of the canal (a . (s).).j 


6 oc (s > 

s 2 0 ci (s) 


6 ec< s > 

<* ci (s) 


-1 


(3,2) 


S 2 + (J)B + (|) 


Available evidence indicates that the system is over- 
damped and ~~ << ~. Therefore equation 3.2 can be approxi- 
mated by 


G ec<*> 



(s) 


-1 


K 


V, 


Cs + S) (s -f £) 


V 


(3.3) 


The short time constant, t « can be calculated from 

s v 

hydrodynamic considerations if the Wavier Stokes equations 
for the canal/cupula system can be solved. Steer (Ref. 65 ) 
solved these equations for a somewhat simplified situation 
and concluded that to first order the short time con- 
stant should be proportional to the endolymph density, the 
square ox the canal's minor radius and inversely proportional 
to the endolymph viscosity. Using the results of Igarashi 
. 38 } for the toroidal radius yields an estimate of 

0,005 seconds for the short time constant in man. 
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The long time constant, has been estimated at 

approximately 10 seconds using subjective responses (Kef, 68 ) 
and approximately 16 seconds using nystagmus records (Ref ., 33 ) 
following step changes in angular velocity. Calculations of 
^ using the audiogyral illusion (Mayne, Ref, 49 ) yield 
values from 8 to 11 seconds. The difference between these 
estimates may be presumed to be due to adaptative processes 
which are more active in the subjective pathways than in those 
associated with nystagmus (Ref, 78 ) , The adaptation dynamics 
will be discussed in the next section# but the important 
point to note here is that neither the subjective reports 
nor nystagmus are merely a consequence of the mechanical 
movement of the cupula described by the torsion pendulum 
model. Therefore any estimate of the long time constant for 
the torsion pendulum model which depends on subjective 
responses must also include the possible effect of neural 
processing. If the presumption were made that no neural 
processing takes place in the vestibular-oeular pathway and 
thus that vestibular nystagmus correctly reflects cupular 
motion, then we would set the long time constant at approx- 
imately 16 seconds. In fact, nystagmus slow phase velocity 
records taken from subjects exposed to steps in angular 
acceleration do seem to show weak adaptation (Young and Oman 
Ref. 62, Malcolm and Jones Ref. 48 ) and thus 16 seconds 
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might be considered a lower bound. Schmid, Stefanelli and 
Mire (Ref. 62 ) calculated the value of by fitting a 
model for the vestibular-ocular dynamics which included an 
adaptation term of the form T A 3 / ( \ s + 1) to nystagmus records 
and found the best fit when ^ = 61.1 seconds and x^ = 18.2 
seconds. Based on this work and its agreement with pre- 
viously cited works, we have chosen a value of. x^ = 18 seconds 
as a good estimate. 

Up to this point, the description of cupular motion has 
been purposely kept vague. It is clear that if the endolymph 
moves within a rigid canal, is incompressible, and no leakage 
occurs around the cupula, then the cupula must move in such 
a way as to sweep out a volume equal to the net volumetric 
displacement of the endolymph. The classical description 
of this movement is that of a swinging motion in which the 
cupula slides freely against the ampula wall and bends near 
its base at the crista. This description was supported by 
Dohlman’s experiments in. which the cupula stained with china 
ink, was observed while pressure was applied unilaterally 


to the fluid and then released (Ref. 23 ) 0 While this 
procedure might indicate that such motion of the cupula is 
possible under application of the pressures employed in this 
experiment, it can not be inferred that such motions occur 
during normal physiologic movement of the head. Steady state 
pressure differences across the cupula have been estimated 
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by Oman and Young (Kef. 59 ) to range from 

1,25 x 10~ 4 dyne/cm 2 at threshold ( 0 1°/ sec 2 ) to approximately 

-*2 2 2 
3.8 x 10 dyne/cm for steps of 3QVsec « Using n short 

time constant of 0.005 seconds and a long time constant of 

20 seconds and presuming a rigid body rotation of the cupula 

about the crista, they calculated a steady state deflection 

2 

of 0.025 degree for a sustained stimulus of 30°/sec . It 
is clear from these calculations, even if they are only 
correct within an order of magnitude that the cupula motion 
observed by Dohlman must have resulted from distinctly non- 
physiological pressures. Oman and Young concluded on the 
basis of these results that the cupula might move angularly B 
iineraly, or both. It should be noted that a linear move- 
ment of the cupula would give superior sensitivity since 
for a given displacement, it would be more effective in 
sending the. sensory hair cells. 

In summary we can conclude that the displacement or 
bending of the sensory hair cells, which is the effective 
agent for iliiciting a change in afferent firing rate, should 
be related to angular acceleration of the head as follows 


Hair cell deflection « 


-a (s) , 


where *r s - 0.005 sec 

x T ~ 18 sec 


(3.4) 


ana 



indicates the inverse Laplace Transform 


Operator 
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A proportional relationship is sufficient for our 
purposes at this point, since the overall gain from head 
acceleration to firing rate can best be estimated from 
records of afferent firing which will be discussed in the 
next section, 

3,1.2 Afferent Processing and Random Signal Variations 

The most desirable data upon which a model of afferent 
vestibular responses in humans could be based would, of 
course, be in vivo recordings in the canals' afferent nerve 
in humans. Since man is not suitable for experimental surgery, 
such data is not and may never be available. There are two 
other sources of data which can be used to make reasonable 
estimates of afferent processes in man. The first of these 
consists of psychophysical data taken from human subjects , 
which of course will include whatever dynamics are present 
in the afferent processing. The second source is recordings 
of peripheral afferents in animals. The sum of data from 
these two sources is not sufficient to argue conclusively 
that a particular dynamic effect is peripheral in man, but 
if such an effect is seen in human subjective responses and 
is also present in the afferent recordings of animals which 
are phy logenically similar to man, then such a conclusion 
seems reasonable , 
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An analysis of vestibular nystagmus and reports of sub- 
jective perception of rotation indicate the need for dynamics 
in addition to the torsion pendulum model to account for rate 
sensitivity and adaptation. The need for rate sensitivity 
shows up principally in cases when there are abrupt changes 
in the rate of rotation. Nashner (Ref ..54, 55) in studying 
human postural control, found it necessary to include a small 
lead term (.017s + 1) in the semicircular canal dynamics to 
predict the response times of subjects exposed to large im- 
pulsive stimuli. A behavior consistent with such a rate 
sensitivity was seen by Benson (Ref. 6 ) in analyzing nystag- 

mus records for sinusoidal stimuli between 0.01 Hz and 5 Hz. 

A consistent increase in amplitude ratio for vestibular 
nystagmus was unexpectedly observed starting at about 0.5 Hz. 
An increase in the amplitude ratio of 3 db is seen at approxi- 
mately 2.6 Hz which would imply a lead term of the form 
(0.06s + 1) . 

The phenomenon of adaptation is much more clearly evident 
in subjective responses than is rate sensitivity. Adaptation 
can be thought of as a fatiguing of sensation which occurs in 
addition to that which arises due to the long time constant 
of the torsion pendulum model. As an example, the torsion 
pendulum model predicts a steady state sensation of constant 
velocity in response to a step in acceleration while sub- 
jective data (Ref. 14,35 ^ indicates a gradual decline in the 
sensation of velocity. For steps in angular velocity, sub- 
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jective data and nystagmus data {Ref. 2 ) indicate not only 
a diminishing of response to zero (consistent with the torsion 
pendulum model) but also a reversing of the response. Finally, 
as noted in the previous section attempts to fit the torsion 
pendulum model to the responses from impulsive velocity 
changes yielded different long time constants for subjective 
and nystagmus data. Young and Oman (Ref. 82 ) were able to 
account for this behavior by adding an adaptation operator 
of the form 


30s 

30s + 1 

to the subjective pathway and 

120s 

120s + 1 

to the nystagmus pathway. The difference, in adaptation time 
accounts for the discrepancy in estimation of the long time 
constant. 

Combining the terms proposed for rate sensitivity and 
adaptation, we conclude that in addition to the torsion 
pendulum model, we should have dynamics of the form 


T 

T 


R 

A 


y ( y + 

(t a s +1) 

W G17 seconds (lead time constant) 


(3.5) 


30 seconds (subjective adaptation time constant) 


sc 
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To decide if it is reasonable to ascribe these dynamics 
to the peripheral afferents, we must resort to data taken from 
the first order afferents in animals. Lowenstein and Sand 
(Ref. 43 , 44 ) and Groen, Lowenstein and Vendrick (Ref* 32 } 
made recordings in the vestibular afferent of the thornback 
ray, Grocn et. al. made recordings in the isolated end organ 
to preclude the possible effects of efferent innervation* 

1 

These experiments confirmed the fundamental features of the 
torsion pendulum model, but make no comment regarding addi- 
tional rate sensitivity or adaptation. The most thorough 
study of the vestibular afferents in a mammal were conducted 
by Goldberg and Fernandez (Ref. 27 ) using the squirrel monkey, 
in these experiments a thorough evaluation of afferent responses 
to constant angular accelerations and sinusoidal stimuli was 
made. The firm conclusion reached was that significant rate 
sensitivity and adaptation was present in a large percentage 
of the cells studied. After eliminating the dynamics which 
can be attributed to the mechanical operation of the endolymph 
cupula system*. Goldberg and Fernandez found afferent dynamics 
of exactly the same form as those given in 3.5. They found 
that t r ranged from .013 to .094 seconds with a mean value of 
.049 seconds. ranged from about 30 seconds to infinity 

(no adaptation) with a typical value being 80 seconds. 

It is clear from these results that the attributes of 
both rate sensitivity and adaptation are present in the peri- 
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pheral vestibular neurons of the squirrel monkey. Since these 
effects, with roughly the same time constants,- are observed 
in human subjective responses we will also consider them to be 
present in the afferent processes of the human vestibular 
system. In addition to the afferent response which depends 
upon the rotational stimulus there is a spontaneous afferent 
discharge and a noise component which is essentially indepen- 
dent of the stimulus. Adding these terms we arrive at the 
model for afferent firing rate shown in Figure 3.1. A very 
conservative figure was chosen for the rate time constant x R 
since its existence in subjective responses is very difficult 
to detect and because it is possible that the most rate sensi- 
tive cells would be used mainly for eye stabilization and con- 
sequently only show up in nystagmus records. 

The constants H and F (see fig. 3.1) have not been deter- 
mined separately but the product HF can be calculated based on 
the magnitude of the afferent response to a controlled stim- 
ulus. Afferent data from Ref. 27 indicates tnat a typical 

initial response to a step in velocity of 1 deg/sec (approxi- 
mately .0175 rad/sec) is about .55 impulses/sec. Substituting 
the values for the time constants into the model and using 
the initial value theorem we obtain a value for HF of -6303. 

The sign of HF can be considered arbitrary as long as the 
processing centers in the brain interpret the sign correctly. 
The value for the spontaneous discharge (SFR) is of little 
consequence for our purposes since it is presumed that the 
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Figure 3.1 Afferent Model of Semicircular Canals 
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higher centers only process differences from the resting 
level, but for completeness we will assign a value of 90 ips 
which is typical of the cells seen by Goldberg and Fernandez. 

<a> 

Lastly, we must address the problem of specifying the 
statistics of the afferent noise, n. The presumption is made^ 
that n is a stationary, gaussian process with zero mean. No 
information is available concerning the autocorrelation of n, 
but the variance of n for different afferent cells has been 
calculated. Goldberg and Fernandez show a histogram for the 
coefficient of variation <CV) for 142 different cells. The 
CV for a particular unit is defined as 


CV = E - ■ ^ i ^ ^ (3.6) 

E{AT) 

* ' 

where AT - time in msec between impulses 

and E{x} denotes the expected value of x. 

CVs varied from about 0.03 to as high as 0.64. The distri- 
bution of CVs showed a sharp peak around CV « 0.06 with two 
thirds of the units falling below CV = 0.25. If t'he higher 
centers were capable of distinguishing between regular and 
irregular units then it would be reasonable to assume that a 
greater weighting would be placed on units with regular dis- 
charge patterns. For this reason, a value of 0.06 will be* 
used as the coefficient of variation typical of the most regular 
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cells found 1/3 of the total population) . Using this value 
we can calculate what one standard deviation in firing rate 
would be for a typical cell. 

E 1/2 (n 2 <t)] = 90 ” = 5,1 ips (3 “ 7) 

Gacek (Ref. 26 ) estimates that there are approximately 
^.2, 000 afferent fibers in the vestibular nerve of the cat. 

Since this would include the otoliths and all three canals 
a figure for one crista of 2400 would be reasonable. The 
equivalent one channel representation of a 2400 channel system 
each with independent additive noise of magnitude o, would be 
one channel with 

E 1/2 [n 2 (t)] = 0//23OO (3.8) 

This reduction in o must be tempered by the following consid- 
erations 

1) We chose CV = .06 which was representative of the most 
regular 1/3 of the total cell population and thus the 
value of 2400 should be reduced to approximately 800. 

2) We presumed in the above analysis that the noise on 
each channel was independent of the noise on the other 
channels. If the noise were exactly the same on each 
channel then there would be no reduction in effective 
noise at all. What the actual correlation might be 

is unknown but it is not unreasonable to assume that 
some correlation exists (especially if the noise were 
related to random movements of the cupula) ■ 
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and finally 

3) We cannot assume that the higher centers are perfect 
in their ability to weed out irregular cells or in 
their ability to combine the resulting regular cells 
in such a way as to minimize the effective noise level. 

In chapter V the noise level necessary to yield 75% 
correct detection for the case of experimentally determined 
threshold stimuli is calculated based on a near optimal model 

i 

of the detection capabilities of the higher centers. This 
results in value of equal to .223 ips which is 

roughly equivalent to 520 independent channels each with a 
noise standard deviation of 5.1 ips. In light of the consid- 
erations listed above this seems to be a reasonable noise 
reduction capability. 

The results of this chapter, to this point, can be sum- 
marized by the following model of afferent firing in response 
to a rotation stimulus: 

Afferent 
Firing = 

Hate 
(ips) 

E 1 / 2 in 2 (t) ] « .223 ips 

u>(s) - rotational rate normal to the plane of the 
canal (rad/sec) . 

. SFR =# spontaneous firing rate (ips) 


f (57.3) 300s 3(. Qls + l),_t (5) + sfr 
.( 188+1) (.005s+l) (30s+iyj 
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3,1.3 Optimal Processing and Model Predictions 

How that a linear model is available for the afferent 
firing rate*, it is relatively easy to develop a model for the 
processing done by the higher centers » An optimal estimate 
can be formulated for any linear combination of the internal 
states of the sensory dynamics if the processor receives 
periodic measurements of the afferent firing rate and has a 
knowledge of the sensory dynamics and the statistical char- 
acteristics of both the measurement noise and the input pro- 
cess. So far we have specified models for both the linear 
dynamics and for the measurement noise but we have yet to 
model the processors a priori knowledge of the stimulus. 

Once this is done the optimal processor can be formulated and 
the entire system including cupula dynamics*, afferent dynamics 
and processing dynamics can be tested and its predictions 
checked against subjective responses. 

It is reasonable to postulate that in most situations 
where a person is exposed to passive motion he has some esti- 
mate of the magnitude of the motion which he expects to ex- 
perience and to a lesser degree an idea of the motions fre- 
quency context. One model for the subject's a priori information 
which incorporates both of these notions quite simply , con- 
sists of a first order filter driven by white noise. This 
is equivalent to modelling the stimulus as an exponentially 
correlated process. The only parameters to be specified 
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in such a model are the filter’s cut off frequency and the 
magnitude of the white noise. Figure 3.2 illustrates the 
higher centers' internal model of the input statistics which, 
together with a model of the sensory mechanism, describes 
the afferent firing which the processor samples periodically. 
One way to communicate to the mathematics of our model that 
very little is known about the frequency content of the stim- 
ulus is to make small and thus render the bandwidth of 
the input spectrum very large. With- this in znind tj. is set 
to a value less than or equal to one second. Q(t) is typi- 
cally set to a constant such that the expected standard de- 
viation of the input process given by Q//2 is essentially 
of the correct magnitude for the actual stimulus being tested. 
Setting Q to a value which correctly represents the stimulus 
magnitude can be justified both on the grounds that one usu- 
ally has a reasonably correct estimate of the magnitude of 
incipient motions, and on the grounds that an approximately 


W(S). 

1 

to (s) 

1 

tj.S + 1 



£{W (t^) W(t 2 ) > = 


E[ft)(t)3 = 0 

If Q(t) ss q then 
£ [a) (t 2 ) ] - 

-2 ~i t l" t 2^ 


2r 


h 


Figure 3.2 Internal Model of the Stimulus Process 
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correct value for Q could be inferred on line based upon the 
level of afferent firing. 

Formulation of the optimal estimate is most conveniently 
presented if the mathematical notation is changed from the 
frequency domain (Laplace transform notation) to the time 
domain (state vector notation) , The processors internal 
moa el of the processes which give rise to the afferent firing 
from the semicircular canals can then be written as follows ; 


where 


and 


x(t) = A X ( t) + B' W{t) 

(3.10) 

y (t) = C x(t) + SFR + n(t) 

x(t) is a state vector which represents the state 

of the canal-stimulus system at time t (4 dim.) 
y(t) is the afferent firing rate at time t (scalar) 
W(t) is the white input process shown in figure 3.2 
SFR is the spontaneous firing rate (90 ips) 
n(t) is the measurement noise. 


The choice of A, B and C used to represent a given linear 
system is not unique. Figure 3.3 illustrates the state space 
realization used here. The standard controllable realization 
(Ref. s ) is used to model the mechanical dynamics of the 
cupula and the dynamics of the h^ir cells. x^(t), x 2 (t) and 
x 3 (t) represent the state of the sensor at time t. x.^(t) is 
xl.o stimulus angular velocity in rad/sec. 




To simulate the process of perception of rotation on a 
digital computer requires first, the simulation of the sen- 
sor's response y(t) to a given input process cu (t) and 
then a simulation of the processing by the higher .centers of 
y (t) • After balancing the need for a simulation with good 

. 'Cr 

:>;>nsory fidelity and the need for reasonable computational 
efficiency it was decided to update the state of the sensor 
every tenth of a second and to update the central processor ' s 

A 

estimate of the rotational rate, {t.) ? every second* The cen- 
tral processor is now faced with the problem of estimating 
m (t) ~ x^(t) given the measurement history of the afferent 
firing rate y(t ) r y(t-l), y (t-2) • . . .y (t-n) • . * » The minimum 
mean squared error (mnse) estimate of x(t) is given by the 
following sequential filter (called a Kalman Filter, Ref. 41 
39 , 76 ) . 


K( t n )[(y(t n )-SFR) ~ SVVV-l^'W] (3.11) 


where £ (t. t , ) is the state transition matrix 

A n, n~ r 

for the system given in 3*10. 
and K(t^.) are the Kalman gains at time t n * 

Since the sensory dynamics are time invariant t n „ i ) 

can be expressed as 61 (t -t , ) and can be calculated as 

J -A n n-1 


follows 
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I A (r) -^{(Is-Ar 1 } (3.12) 

where A is defined by equation 3.10 and figure 3.3 

A 

The optimal estimate x{t n ) differs from the true state 
x(t ) by an error £(t n ) which is zero mean and has a covari- 
ance given by: 


P(t n ) = E[£(t n ) £ T (t n >] = (I - K(t n )C)P' (t n ) (3.13) 

where ]? 1 It ) is the covariance of the processors know- 
ledge of xft^) given all past measurements of the afferent 
firing y(t n-1 ), y(t n- _ 2 ).,. but not y(t n ). P’ (t R ) is given by 

£' <V = Ia'VVi^^-iL^VVi' ■ 


^ t n" t n-l ) 

J A (TJ B d"(r)B T f^(r)clr 


(3.14) 


The Kalman gains K(t n ) are calculated from P'(t ft ) as 

follows: 

K(t n ) = £' (t n )C j'c£'(t h )C T + E[n 2 (t n ))j (3.15) 


Although a full exploration of equations 3,11 - 3.15 
cannot be given here, some motivation for the form of equa- 
tions 3.11, 3.14 and 3.15 can be given. 

Equation 3.11 shows that the minimum mean squared error 
(mmse) estimate of x(t n ) is made up of two terms. The first 

A 

term £(t R )x(t n _-, ) merely propagates the optimal estimate 

at time t .. forward in time and represents the state which 

n-JL 

A 

would exist at time t_ if the estimate x(t n _^) were errorless 

and the system had no inputs during the interval (t , t ) . 

i * — x • n 



56 


Alternatively this term can be thought of as the best estimate^ 
of x(t n ) based on the measurements y(t n _^), y (t^^) 11 0 0 0 
The second term represents the difference between the ex- 
pected afferent signal ^E[(y(t n > ~ SFR) /y ( fc n ^) t ■> • » • I = 

CO(t ,t , )x(t n)l based on the old measurements and the new 
-“A n n-J. n-r j 

measurement of the afferent signal (y(t n ) - SFR) * This 
second term therefore summarizes the relevant new informa- 
tion from the latest measurement. Note that the afferent 
signal is always measured from the spontaneous firing rate 
(SFR) since SFR is independent of the state x (t n ) « With 
these interpretations of the terms in equation 3,11 the 
Kalman gains K(t R ) can be interpreted as the weighting factors 
which indicate the relative importance or usefulness of the 
new information as compared to the old information in estimat- 
ing the state vector at time t . 

P s (t n ) given by 3.14 represents the error covariance of 

an estimate of x(t n ) based only upon measurements taken before 

T 

time t ♦ The term 4> (t ~ t ~ ) P ( t„ -a ) (t -t i > represents 
the covariance of the estimate error £’ (t n ) * x(t n ) - 

- (t n, t n-l^ due to the error = ~ 

x ( t j ^ _ i ) propagated forward in time. The integral term rep- 
resents the covariance of §*(t n ) arising from the unknown 

stimulus during the interval (t . t. ) . 

n— x f n 

Roughly speaking equation 3.15 can be looked upon as the 

/ 

ratio of the variance of the old information divided by the 
sum of the variance of the old information plus the variance 
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of the new information. Note that the leading term is not 

quadratic in c since the term which Ktt^) multiplies in 

equation 3,11 is already linear in C, Thus K(.t ) increases 
2 

if E[n (t n ) 3 decreases and therefore weighs more heavily a 

measurement with a smaller noise component. Alternatively 

K ( t ) decreases if P’ (t ) decreases since this indicates that 
n n 

the accumulated old information is relatively more useful 
than the new information gained from y(t n ). 

Implementation of equations 3.11-3.15 is relatively 
simple once it is recognized that most of the expressions can 
be calculated in advance. For the simulations carried out 
here measurements of the afferent firing rate are made avail- 
able once every second which implies that t = t n - = 1 

second. Equation 3.12 is used to calculate $ ( t ». 1) and the 
result ( a 4 X 4 matrix of constants) is stored for future 
calculations. The integral term in equation 3.14 (call it 
Pj(t n )) is also a 4 X 4 matrix of constants which can be calcu- 
lated if Q(t) is known (see Figure 3.2). 

To start the simulation one must have an initial state 

-a 

estimate x(t Q ) and an associated error covariance for that 
estimate, P(t Q ) . Between time t Q and the first measurement 
of afferent firing at time t^ the processor computer P' (t^) 
from $ A (i), P(t Q ) and P^tt^) and then the Kalman gains Kft^) 

from P l (t^) , C and the variance of the measurement noise 

2 

E[n (t^)]. Since we have only one measurement the inversion 
implies simple division. Once the measurement y(t^) is avail- 
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able equation 3,11 can be used to calculate the new estimate 

A 

^(t^) and equation 3,13 can be used to calculate the associated 
error covariance P (t^) , While waiting for the measurement 
y(t 2 ) the processor repeats the above steps calculating P # (t ? ) 
and K(t 2 ) * etc. 

In a real time application in which the measurement noise 
variance K [n (t^) ] and input power Q(t) are known in advance 
it is possible to precalculate P* (t ) Q K(t ) and P(t ) for 
all future times t n , since they are independent of the measure- 
ment y(t n >. For systems which are asymptotically stable and 
for which the measurement noise n(t^) and input power Q { t) = Q 
are time invariant P* (t R ) , K (t^) and P (t n ) approach constant 
values (denoted and P^) . In such a case equation 3,11 

becomes time invariant (upon substituting K for K (t ) ) and 

A A 

represents the state space version of the Weiner Filter 
(Kef ,60 / 73 ) , 

Appendix I lists the calculated values for the transition 
matrices, Kalman gains, etc., which pertain to the simulation 
of perception using afferent information from the semicircular 
canal system. 

To test this model of suprathreshold vestibular percep- 

2 

tion a step in angular acceleration of l,5°/sec was simulated. 
The stimulus was on for 120 seconds and then off for 120 sec- 
onds leaving a constant rotational velocity of 180 degrees/ 
second. The response of the model (shown in Figure 3,4} peaks 
approximately 27 seconds after the onset of the stimulus and 
decreases to less than 10% of the peak response after two 
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.0262 t 
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0 < t < 120 
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Time (seconds) 


Figure 3 . 



4 Predicted Subjective Response to 1.5 Degree/Second 
Acceleration Step 
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minutes . When the acceleration ceases the predicted percep- 
tion quickly changes sign with a secondary peak occuring 28 
seconds after the stinulus is removed „ 

Clark and Stewart (Ref. 14 ) and Guedry and Lauver (Ref, 

35 ) have conducted experiments to determine the subjective 

<> 

' 2 

response for acceleration steps of 1.5°/sec . Although the 
data from these two experiments agree in a qualitative way, 
there are significant quantitative differences between them. 
Clark and Stewart show a peak response which occurs at approxi- 
mately 35 seconds with essentially no response at 120 seconds 
after onset. After the acceleration is removed (at 120 seconds) 
a reversed response with essentially the same time course but 
diminished in magnitude is reported. It should be noted that 
a linear model will not predict a reduced secondary response 
if the original response has stabilized at zero. 

Data from Guedry and Lauver show a peak response at about 
25 seconds and an adaptation time constant of approximately 

30 seconds (see ref. 35 ) 0 One interesting aspect of tnis 

/ 

data is that it indicates an initial rate of change of sub- 
jective velocity equal to approximately 2.9 (degrees/second)/ 
second. This is almost twice the true rate of change and 
would imply that the initial response to a step in angular 
velocity would also be twice the true rate. 

The response predicted by the model reaches an initial 
peax at 27 seconds (consistent with Guedry and Lauver ) , ap- 
proaches a negligible response at 120 seconds (consistent 
with both sets of data) , has a roughly symmetrical response 
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after the step in acceleration is removed (dictated by its 
linearity) and has a gain which yields an initial perception 
of angular velocity consistent with the true stimulus magni- 
tude. 

In addition to the model's ability to predict subjective 
sensation for rotations about an earth vertical axis, it has 
a degree of flexibility not found in earlier models. Since 
all of the dynamic characteristics observed in human subjec- 
tive responses are also observed in the afferent recordings 
made in animals the model assigned these characteristics to 
the human afferent. Consequently it was necessary to postu- 
late a high bandwidth stimulus process so that the optimal 
filter would not contribute significantly to the overall dy- 
namics of tne subjective response. If it were discovered that 
the afferent dynamics in man were significantly different 
than tne dynamics observed in man's subjective responses then 
the optimal estimator model for the higher centers might con- 
tribute to an explanation for the discrepancy. One example 
of sensory processing which contributes significantly to the 
overall dynamics of perception is the processing of afferent 
signals from the otoliths which is described in the remaining 
sections of this chapter. 
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3 0 2 Otolith System 

Specification of the afferent dynamics for the otolith 
sensors differs in one fundamental respect from that of the 
semicircular canal system. In the case of the semicircular 
canals, definative data was available for the response of 
mammilian afferents to dynamic stimulation while no similar 
data is presently available for the otoliths. In this sec- 
tion the afferent response of the otolith sensors will be 
inferred from the . following information 

(1) The general mechanical structure of the otoliths 

(2) The data available on the static response of the 
otoliths to a constant sheer force 

(3) The limited afferent data available for dynamic 
stimulation plus some qualitative comments re- 
garding the nature of the afferent response seen 
in the squirrel monkey by Dr. Goldberg (personal 
communication) 

and (4) Tne known subjective response to accelerations from 
human subjects 

coupled with the presumption that the higher centers process 
tne afferent data optimally to yield a perception of specific 
force. Specific force is defined to be £ - a where is the 
local gravitational force vector/ and a is the acceleration 
of the body with respect to inertial space. Stated differ- 
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ently, specific force is the gravitoinertial reaction force 
per unit mass. 

3.2,1 Division of Afferent Response ahd Higher Order Pro - 
cessin g 

The basic mechanical structure of the otolith sensor can 
be modeled as a mass (the otolithic membrane) immersed in a 
fluid (endolymph) and restrained by the spring-like action of 
its mechanical restraints. Movement of the Otoconia which 
is presumed to be proportional to hair cell displacement can 
thus be described as follows: 

Mx(t) = -Kx (t) - Vx(t) + M SF(t) (3.16) 

where x(t) is the displacement of the otoconia from 
their resting position 

M is the mass excess of the otpconia over an 
equal volume of endolymph 
K is the effective spring constant of the sys- 
tem per unit displacement of the otoconia 
V is the coefficient of mechanical and viscous 
losses due to movement of the otoconia 
and SF is the specific force acting on the otoconia 
parallel to the plane of the macula. 

Using Laplace transform notation we can solve for the 
transfer function relating otolith displacement to specific 
force : 


X (s) 
SF(s) 
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(3.17) 


2 , V 
S + M S 


K 

M 


An order of magnitude estimate can be made for the nat- 
ural frequency w n = based upon estimates of the displace- 
ment of the otolith organ under the influence of a lg shear 
force (the following procedure for estimating is due to 
Young Ref. 03 and through personal communication with C. Oman). 
Applying the final value theorem to equation 3*17 with 

SF(s) = — we find that the steady state displacement of the 
s 

3 2 . 

otolith organ for a lg (10 cm/sec ) stimulus is 

x<«) = £ 10 3 cm (3.18) 

ft 

Two rough estimates can be made for this displacement. The 
first results from assuming that the sensitivity of otolith 
hair cells is approximately the same as that for the semicir- 
cular canals and thus the displacement of the otoconia must 

be approximately the same at threshold as the threshold dis- 

placement calculated for the cupula. Oman (Ref„ 58 ) calcu- 
lated the threshold displacement of the midpoint of the cupula 
to be approximately 10~ 6 cm. If this were the displacement of 
the otoconia at threshold (approximately .005 g) then a lg 

stimulus should produce a displacement of x(») = 2 x 10 4 cm. 

-7 

Usrng this figure in equation 3.18 yields a value of 2 x 10 
M Pk 1 

for ana thus w n = f ^ =* 2240 rad/sec. The second estimate 
for x(«) is given by deVries (Ref. 20 ) based upon X-rays of 
the ruff sacculus. He found a steady state displacement of 
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~ 3 

x(°°) ~ 7 x 10 cm for a Ig stimulus. This value for x(°°) 
results in an estimate of W n = 374 rad/sec. 

Although these estimates differ by a factor of six they 
both confirm that the natural frequency of the sensor should 
be at least an order of magnitude above the highest frequency 
normally associated with vestibular stimuli. With this estab- 
lished we can now demonstrate that the transfer function in 
equation 3.17 can be represented as either a pure gain or at 
most a simple lag over the frequency range of interest, 

u)e (° • 

10 

Figures 3. 5A and B show the Bode amplitude ratio plots 
and phase angle plots respectively for equation 3.17 for dif- 
ferent damping rations £ = V/(2v^lK) £ 1.0. If £ £ 1 then it 
is clear that for frequencies less than w n /lC the system can 
be approximated by a gain. If £ > 1 then the denominator of 

equation 3.17 has two real roots and s 2 which satisfy s-^s 2 
2 

^ n which implies that one root must have a value greater than 

u> and one root a value less than to . It is only the one root 

with a value less than uj n which could significantly affect the 

frequency response of the system for uj < w n /10. Thus for 

a) £ (0/° n) , which should include all. of the normal physiologic 
10 ' 

motions, we can approximate the mechanical dynamics of the 
otolith organ by: 

x (s) K A 

SF (s ) 3 + A (3.19) 




Figure 3.5 Amplitude and Phase Plots for 
Second Order Systems 
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where A = ~~ + j(^j " if S > 1 

M 

and A = 00 if £ < 1 

The transfer dynamics from otolith displacement to 
afferent firing must now be considered, Vidal et al (Ref. 70 ) 
in experiments with cats found a unidirectional rate sensi- 
tivity of the form 

Firing rate - G(0(t) + K0(t)3 (3.20) 

where K = 0 0(t) <_ 0 

K > 0 §<t) > 0 

and 0(t) is the tilt angle 

This response nonlinearity can be removed if it is assumed 
that perception is a function of the difference in firing rates 
from cells with opposing gains. Such an assumption is reason- 
able based upon the work of Flurr and Mellstrbm (Ref. 25 ) and 
leads to the following response to tilts: 

Firing rate = G(20(t) + K0(t)) (3.21) 

regardless of the sign of 0(t). We have thus retained the 
effect of rate sensitivity while eliminating the nonlinearity. 
Vidal postulated the model given in 3.20 based upon a stimulus 
whose frequency content was confined to w <,i Hz and thus 
might not be expected to indicate the mechanical dynamics pos- 
tulated in equation 3.19. In the static experiments conducted 
by Vidal no adaptation was reported even though recordings 
were made for as much as 3 minutes at a given orientation. 



Fernandez and Goldberg (Ref, 24 ) indicated that the 
majority of cells from the otolith organs in the squirrel mon- 
key achieved an essentially constant rate of firing within 30 
seconds of a position change. Thus although some cells may 
exhibit significant adaptation a large percentage of the pop- 
ulation does not. 

While it is recognized that the otolith afferents which 
exhibit adaptation may contribute to the transient response 
(in exactly the same manner as nonadapting ©ells) they do not 
govern the long term sensations due to the otoliths which do 
not show significant adaptation to static tilt angles. There- 
fore we postulate the following form as a model for the re- 
sponse of a non adapting otolith afferent to the component 
of specific force in the plane of the macula % 


(^Afferent *) 

) Firing Rate J (s) 


[ — .. .. ± a - ]sf( 8 ) + SFR + n (s) (3.22 

s + A ' s 


where SFR denotes the resting discharge 

n(s) is the measurement noise process 
and A,B,C are shown in figure 3,6 
If, as in the case of the semicircular canal model, we 
postulate that the higher centers’ a priori information con- 
cerning the input (in this case SF) is given by white noise 
through a first order filter and that the higher centers pro- 
cess the afferent information optimally to estimate SF then 
we can view the entire sensory process for the otoliths as 
shown in figure 3.6, 
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Figure 3.6 Model of Otolith Perception 
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The steady state optimal processor H° (s) can now be de- 
termined as a function of W f & s C and N . Since the 

resting discharge SFR is independent of the input and is known 
by the higher centers it is presumed that it is subtracted 
from the afferent signal and plays no further part in the 


processing., 

To find H° (s) we must solve the associated Wiener Hopf 
equation. Since H° (s) must be causal the associated Wiener 
Hopf equation is most easily solved by the method of spectral 
factorization (Ref. 60 ) . This yields the following solution 
for h 0 is) ; 


H° (s) = K' 


S+A 

(S+F) (S+G) 


Where F , G and K° must satisfy 

1 


f 2 + g 2 - bV 


2 2 

V N 


+ A' 


2-2 (B+C) 2 A 2 W 2 , A 2 

F G » r — 5 + 7Z 

y N \ 


W 2 [(B+C) A + ~ j 

j o __ 1 

^ 2 2 1 1 

T ■ TT(F+^> (G+~~) 

1 T I 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


Now that the form of H° (s) is determined we can cascade 
it with the mechanical and afferent dynamics (given in 3.22) 
to yield a prediction for the form of the overall perceptual 
dynamics associated with the otoliths. 
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{Mean Perception of SF) _ Bs + (B+C)A yo s+A 

SF ( s ) * s+A (s+F) (s+G) 


BK° 


(B + ^±|)A) 
Ts+Fj ( s+G) 


(3.27) 


Equation 3.27 represents the framework within which we 
must be able to predict subjective perception which involves 
the otolith organs. Young and Meiry (Ref. 81 ) proposed the 
following model for the perception of acceleration*. 


{Pe r ceived Lateral A cc eleration} (s) _ 1-. 5(s+ .076) , ~ ?R v 

Ti’rue Lateral Acceleration } I (s) “ (s+.19) (s+1 . 5) ' 

It is clear that equations 3.27 and 3.28 have identical 
forms but it is not clear that the parameters (W, , A, B, C 

and N) can be chosen in such a way that the ooef f icients in 
the present model will be in essential agreement with those 
in 3.28. 


3.2.2 Ot olith Model Spe cifica t ion and Predictions 


In this section the criteria for specifying the models’ 
parameters will be elaborated. Based on these criteria values 
for the parameters will be chosen and the resultant model 
evaluated. Actually two sets of parameters must be developed', 
one for the utricle and one for the saccule. The only dif- 
ference is that the resting discharge rate, the sensor gain 
and the measurement noise level for the saccule have been 
found in the squirrel monkey (Ref. 24 ) to be approximately 
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half those found for the utricle. Thus we have: 


Utricle 


A 

B 

C 

N 

SFR 


Saccule 

W 

T I 

A 


C/2 
N/2 
SFR/ 2 


( 3 . 29 ) 


Inspection of equations 3.24, 3.25 and 3.26 reveals that 
F and G will remain unchanged by this transformation while 

K saccule wil1 be twice the value of Utricle t0 make up for 
the difference in sensor gain. With this in mind we will pro- 
ceed with the specification for the utricle model and then 
specify the saccular model by employing the above transforma- 
tion . 


Fernandez, Goldberg and Abend (Ref. 24 ) found an average 
steady state change in firing rate from the utricle due to 
a lg step in specific force to b e 45ips« If the final value 
theorem for a lg specific force is applied to equation 3.22 
we find the steady state response to be B+C. Thus, our first 
condition is that 

B+C = 45 (3.30) 

Since our model of subjective sensation (equation 3.27) 
must fit the same data from which the Young and Meiry model 
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(equation 3.28) was constructed, it is reasonable to assume 
that the corresponding coefficients should be approximately 
the same. It should be noted though that any combination of 
differences between these parameters which preserve the essen- 
tial. fit of the model to the data would be perfectly acceptable. 
Setting corresponding coefficients in approximate equality , 
yields the following four conditions* 

(B+C)§ = .076 (3.31) 

F = .19 (3.32) 

G -1.5 (3.33) 

BK° *1.5 (3.34) 

Equations 3.32, 3.33 and 3.34 should be substituted in 
equations 3.24, 3.25 and 3.26 to yield. three nonlinear algebraic 
conditions in W, , A, B, C and N while eliminating F, G and K° 
Once all the parameters of the model are chosen it will be 
possible to calculate the, afferent signal to noise ratio (S/N) 
which the. processor has assumed in arriving at its estimate of 
specific force. If an unreasonable signal to noise ratio must . 
be postulated to satisfy the other conditions that we have set, 
then the hypothesis of an optimal processor would have to be 
seriously questioned. The signal to noise ratio could con- 
ceivably be set in one of two broad regions. Either the pro- 
cessor would be optimally structured to interpret signals near 
threshold when s/N ~ 1 (this choice would be best suited to 
tasks like the control of postural sway) or it could be struc- 
tured to best handle large signals for which s/N is signifi- 
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cantly greater than one. The only case which would indicate 
a serious problem would be a case in which s/N << L There- 
fore, after calculating the signal to noise ratio at the affer- 
ent level in terms of the model parameter we can set the 
following restrictions 


Ej (Al'R—S FR) ] 
E[n 2 ] 


E s/N 85 


w 


B 2 /t j-MB+C) 2 a 
# 

2{A+i 

T I 1 


1 or greater 

(3.35) 


As in the case of the semicircular canals, there is data 
available concerning the noise level on the afferent fibres. 
Fernandez et a.l (Ref. 24 ) plotted the CVs obtained from 47 
otolith units and concluded that while otolith CVs are signifi- 
cantly lower than canal CVs there is no difference between the 
CVs of the utricle and saccule. Since the level of spontaneous 
firing at which the utricle afferents operate ( =88ips) is 
twice that for the saccule ( -‘44) we conclude that using the 
same CV will yield a value for N saccale \tricle ^ 2o In “ 
spection of the CV histogram offered by Fernandez indicates 
that a value of CV ~ .04 would be representative of the most 
regular afferents. Using a tonic firing rate for the utricle 

of 80 we calculate the standard deviation of the firing rate 
1/2 2 

to be E [N ] = 3.39 ips. In the case of the semicircular 
canals we found that the central processor was capable of re- 
ducing the effective measurement noise by a factor of 1//520 
by utilizing the large number of individual afferents avail- 
able from each sensor. If we presume as a first approximation 
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that the central processor for the otoliths is equally capable 
of combining afferent signals to reduce the effective measure- 
ment noise, then we could expect; 


E 1/2 [n 2 ] = N 


3.39 


= .149 ips 


(3.36) 


Finally it is reasonable to assume that the presumed 
bandwidth of the input, which is governed by should be 

of the same order of magnitude as that found necessary for 
the semicircular canals. If there is a difference, then 
1 /t.j- for the otoliths should be slightly lower than that 
for the canals since the otoliths are generally considered 
to be sensitive to lower frequency stimuli than are the 
canals. Thus we should find 1/x^. roughly in the region: 


.5 rad/sec = 


5 rad/sec 


canals 


'otoliths 


canals 


(3.37) 


After substituting the values for F, G and BK° from 
equations 3.32, 3.33 and 3.34 into equations 3.24, 3.25 and 
3.26 we have a set of eight algebraic conditions (3.24, 3.25, 
3.26, 3.30, 3.31, 3.35, 3.36 and 3.37) to be satisfied by 
six parameters (W, A, B, C and N) . The following values 

for the six parameters were found to best fulfill the stated 
conditions (the value of the resting discharge rate is given 
for completeness) . 




76 


Parameter 



Utricle 

.00268 


Saccule 
. 00268 


o 2 


o 2 


B 

C ‘ 


91.1 45=55 

46.1 -23.05 


N 

SFR 


o 147 


88 


.073 


44 


The consequences of using these parameter values are as 
follows x 


(1) B*fC = 45 (fulfills eqn 3*30) 

(2) Equations 3.24 and 3.25 are fulfilled exactly but with 

F = ,133 (see 3,32) 

G = 1.95 (see 3.33) 

(B+C)| = .0988 (see 3. 31) 


Each of these values represents a 30£ change from the approxi- 
mate values specified. Since these values represent the critical 
frequencies of the overall subjective model, their only import- 
ance arises in how they affect the models prediction of the 
system's phase response. Figure 3.7 shows the phase plot for 
Young and Meiry's model (eqn. 3.28) and the model derived here 
(egn, 3.27 with the parameter values specified above). The 
phase predictions of the model are quite good over the fre- 
quency range above .5 rad/sec and only slightly lower than 
the few data points available below .5 rad/sec* 




Figure 3.7 Comparison of Phase Predictions for Otolith Models 
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(3) BK° = *911 (reasonable agreement with. 3,. 34) 

(4) s/N = 1.21 (sec 3*35) 

(5) N * o 147 (sec 3«36) 

and (6) = 1 (sec 3 0 37) 

The resulting afferent dynamics are given by 


Afferent I 
Firing V (s) 
Rate ) 


91,1s + 9 

st o 2 


SF (s ) + 


+ n (s) 


(3.38) 




where E^ 2 [n^(t)J = .147 

Figure 3.8 shows the change in afferent firing due to a 
step in acceleration of Ig. The response consists of an im- 
mediate jump of 91„1 ips which decays with a 5 second time 
constant to a steady state change in firing of 45 ips„ It 
should be stressed that this five second decay could be due 
to either the mechanical dynamics (implying that the mechan- 
ical dynamics are very over damped) or could arise from the 
dynamics of the afferent processes „ If more phase data be- 
comes available on the subjective response to stimuli below 
o 5 Hz and confirms the greater phase lead shown by Young and 
Meiry, the model can be altered to give an identical phase 
prediction. The penalty for this change is twofold. First 
the resulting model for the afferent response shown in figure 
3.0 would jump immediately by 400 ips and then fall to 45 ips 
with a time constant of approximately 1.5 seconds (which at 
first glance seems less reasonable) . Secondly,, the resulting 
prediction for acceleration step thresholds would be signifi- 
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80 

cantly reduced below currently accepted figures. 

In a personal communication with Dr* Goldberg it was 
learned that his current research with the otolith organs of 
the squirrel monkey shows a wide variety of afferent dynamics, 
some of which are consistent with the general prediction of 
Figure 3,8, When further experimental data is available con- 
cerning the dynamic response of otolith afferents an analysis 
similar to the one made for the semicircular canals can be 
undertaken. In the meantime it is interesting to note that 
the approach taken here can yield a model which accounts reason- 
ably well for the available subjective data, the known physio- 
logical structure of the sensor and makes reasonable predic- 
tions concerning the afferent processes and the associated 
central processing. 

Since for time invariant systems the steady state Kalman 
filter is equivalent to the Weiner filter developed above, it 
was decided to use the time domain formulation of the problem 
for the otoliths as was done for the semicircular canals. 
Appendix I gives the appropriate transition matrices and Kalman 
gains for the otolith simulation. Figure 3.9 shows the result- 
ing subjective response of the overall system to a lg step in 
acceleration. While this prediction shows less overshoot than 
the model by Young and Meiry there is no subjective data with 
which to directly compare it. 
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Chapter XV 

QUALITATIVE NATURE OF PROCESSOR FOR DETECTION OF NEAR- 
TIiRESHOLD HORIZONTAL ROTATION 


Numerous studies have been made to determine the magni- 
tude of vestibular thresholds in the three rotational axes 
(with primary emphasis on the yaw axis) . A comprehensive re- 
view of these experiments through 1965 is given by Clark (Ref. 
10 ) . It is interesting to note that numerous definitions 
of threshold are used by experimenters and thus it is diffi- 
cult to compare the experimental results. An even more vexing 
problem is that the results of these experiments cannot be 
generalized to predict the probability of detection as a func- 
tion of time for an arbitrary rotational stimulus. The reason 
that such a prediction cannot be made is that previous inves- 
tigators have not proposed a stimulus-perception model ade- 
quate for arbitrary near-threshola rotations. In this cnapter 
two fundamentally different hypotheses are considered for the 
nature of the threshold mechanism. The consequences of each 
hypothesis are explored and an experimental procedure designed 
and carried out to test their validity. 

4 . 1 Processor Hypotheses 

When a man is rotated, what information is available 





82 

which he can use to estimate the time course of his rotation? 

If the subject’s eyes are open the images of a stationary 
environment moving across the retina provide one source of 
information- For those cases in which the visual environment 
is not inertial ly stable (for example, it may rotate with 
the subject) or in which the subject's eyes are closed, another 
input is needed- Although cues arising from movement of the 
bodily organs or from tactile sensations may be present, the 
most important non-visual source of information concerning 
rotation arises from the vestibular sensor- If a subject is 
denied visual cues and tactile and proprioceptive cues are . 
ignored, then the higher centers are left with the information 
provided by the vestibular afferents in the form of change's in 
neural firing rates. 

Dynamic models for the mechanical movement of the cupula 
and consequent afferent discharge were described in C hapter 
Throe. Since these models were developed as part of a model 
of suprathreshoid perception no consideration was given to 
either a mechanical or neural threshold nonlinearity. The 
existance of thresholds for rotational stimuli is sometimes 
accounted for by a nonlinearity at the sensor which would imply 
that until a certain stimulus level is reached there is no 
stimulus related cnange in afferent firing level. It is this 
general conception of the threshold mechanism (shown in F igure 
4.1) which will be called the "simple threshold model". The 
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Figure 4.2 Signal in Noise Model 



particular form of the nonlinearity is arbitrary except for 
the dead zone . 

The second model will be called the "signal in noise" 
model since, as shown in Figure 4 = 2, it consists of nothing 
more than a measurement noise source between the sensor dynam- 
ics and the processing by the higher centers. The threshold 
phenomenon in this case arises from the masking effect of the 
measurement noise at low signal levels. The magnitude of the 
measurement noise is set so that a correct detection of the 
direction of motion will result seventy-five percent of the 
time . 

Although these two models represent two fundamentally 
different approaches to modelling the phenomenon of vestibular 
thresholds (one including and one excluding an afferent signal 
correlated with subthreshold stimuli) they both predict the 
general variation of response latencies as a function of stimu- 
lus magnitudes. Figure 4.3 illustrates the time history of 
the cupula displacement (assuming a linear mechanical re- 
sponse) for steps in angular velocity and acceleration . For 
steps iii angular acceleration, the displacement builds up to 
a steady state value with a time constant of IB seconds. 
Applying the "simple threshold model" we see that if the 
steady state displacement is less than then the output 

of the nonlinearity will remain zero and the performance of 
the processor in determining the direction of the stimulus 
should remain at 50%. The larger the steady state displace- 
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merit of the cupula, the sooner the output of the nonlinearity 
deviates from zero, which results in shorter latencies as is 
observed in experimental data (this is essentially a descrip- 
tion of the classical cupulogram, Ref. 68 ) ° In the case of 
a velocity step input, detection will generally come quickly 
if the input is above threshold since the peak cupula displace- 
ment follows the stimulus onset with almost no delay. 

The ’’signal in noise" model also predicts similar changes 
in latency. If the signal to noise ratio is taken as a measure 
of the information available to the higher centers, it is clear 
that for the case of an acceleration step input, the informa- 
tion flow peaks at about 20 seconds while in the case of a 
velocity step the best signal to noise ratio occurs immediately 
after the onset of the stimulus. If the stimulus magnitude 
is increased, a sufficient amount of information to make a 
response becomes available earlier in the case of an accelera- 
tion step, while the latency will continue to be short in the 
case of a velocity step. 

4 . 2 Threshold Characteristics of the Hypothesized Models 

Since both models are capable of explaining the existence 
of thresholds and the general trend of response latencies, it is 
necessary to design an experiment in which the results pre- 
dicted by the two models are measurably different so that one 
or both of the models can be rejected. 
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The maximum cupula displacement for the velocity step 

V(t) and the acceleration step A(t) shown in figure 4.3 are 

the same. If this maximum value were slightly greater than 

D MIN 4.1) then the ’’simple threshold model” would 

predict that both V(t) and A(t) would be threshold level 

stimuli. Combining these velocity step and acceleration step 

stimuli at their threshold level results in a stimulus which 

would produce a step in cupula displacement just larger than 

D...... . Therefore if we define 

MIN 

u> (t) = V{t) + A(t)t (4,1) 

then the "simple threshold model" would predict that w (t) 

C 

would also be a threshold level stimulus. 

This result does not follow for processors which deal 
optimally with signals in noise. If the noise- free firing 
rate due to w (t) is compared to either that for the velocity 
step threshold or. the acceleration step threshold it is clear 
that with the same level of additive noise the presence of 
co { t ) is more likely to be detected than either V(t) or A(t) 
alone. If the signal in noise model is correct, co (t) will 
have to be reduced in magnitude if it is to be a threshold 
level stimulus. The questions which naturally arise are 

1. By what factor must u> (t) be reduced? 

c 

and 2. Is there a sufficient difference to detect experi- 
mentally? 

A measure of the information content of a signal s{t) 
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at a given time? when immersed in independent zero mean addi- 
tive gaussian noise, n(t), is the signal to noise ratio. S/N, 
defined as 


5 

N 


E { S 2 (t) } 

E {nht) } 


(4,2) 


where E {°} indicates expected value * 

Let n(t) be a stationary process and define N = E{n (t)} Q 
Since, in the problem at hand, the signal is available over a 
period of time, the integrated signal to noise ratio is the 
appropriate measure or: information content 0 If FR^ft) and 
FR (t) denote the afferent firing rates due to threshold level 

A 

velocity and acceleration step stimuli respectively, then the 
"“signal in noise model requires that 



r T cra A 2 it» dt _ 

J = 7 


or that 


T 

J FRyhtjat = 

6 


FR (t)dt 

A 


(4,4) 


where (Q,T) is the interval given to the subject to detect 

the stimulus (usually 10-15 seconds),, The problem is to find 

a factor, K f such that Kw (t) gives rise to a firing rate, 

c 

call it FR (t) , which has an information content equal to I* 


Since 
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FR (t) = K (FRyCt) + FR A Ct)) (4.5) 

we can calculate its information content a3 follows? 


T T 

f FR c 2 (t)dt = K 2 J (FRytt) +FR A (t)) 2 dt 


* K 


FR^ 2 (t) dt + J FR A 2 (t)dt + 2 J FR^tjFI^ 


FR„ (t) dt) 


~ 2K \ I + 


FR v (t)FR A (t)dt 


(4.6) 


If we set this equal to I we can solve for the required Ks 

1/2 


K = 


/2 l I + J^FRyCtjFR^tJdt 


(4.7) 


Several facts should be noted at this point? 

1. K is independent of the noise level as long as n(t) 

is stationary, but is implicitly a function of the 
sensor dynamics, ' through FRy(t) and FR^ (t) . 

2. For all reasonable choices of 

J FR v (t)FR A (t)dt > 0 (4,8) 

o 

since FR y (t) and FR A (t) share the same sign in the 
interval (0,20) seconds. 

3. Finally, using the Cauchy Schwartz inequality, we 


have 
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T 

/ FRy (t) FR^ (t) dt 

O 


< 



aa 



(t) dt 


1/2 


(A. 9) 


Combining 4.8 and 4.9 we can conclude that 

T 

0 < J FRyttJFR^ttJcit < X 

O 

which, utilizing 4.7 implies that 

1/2 < K < 1//2 


(4.10) 


(4.11) 


If instead of integrated signal to noise the above 

analysis used a first order lag) 1/(ts + 1) ; f or which we will 

C T Lfx}dt 

use the notation where L{x} denotes the response of 

the linear system s/(ts + 1) which is initially at rest to the 

input x j the results would be essentially unchanged. Specifi- 
cally we haves 


P P 

[ LlFR^Jdt = j L{FR A 2 }dt 


(4.12) 


and therefore 


L{FR c 2 }dt = j" L{K 2 (FR V +FR A ) 2 }dt 

2 f T 

= 2K (1+ j L{FR v FR A )dt) (4.13) 


Since. 


r T 

J L {FRyFR^} dt > 0 


and 


J L{FR y FR A }dt < ^| T L{FR v : 


1/2 


}d- 


T 


(4.14) 

1/2 


L{FR a }dt 


= I 


(4.15) 
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we can again conclude (after setting 4.13 equal to I) that 

1 < K < - (4.16) 

2 ft 

In fact the above analysis is valid for any linear 

f T 

operator, l L{-}dt, as long as it satisfies the criteria 
required of an inner product. This result will be useful 
in Chapter Five when we consider models for the detection 
process . 

Figure 4.4 summarizes the above predictions. Every 
point in the figure represents a stimulus made up of a 
velocity step component v and an acceleration step compon- 
eat a. For example the point P shown in the figure represents 
the stimulus 

P(t) « (a + St) deg/sec (4. 17) 

The points labeled V and A represent the velocity and 
acceleration step threshold levels respectively for a given 
subject. The dotted line through the origin and the point 
(A , V) represents the class of stimuli over which the two 
models differ in their prediction of threshold levels most 
dramatically. The threshold prediction of the "simple 
threshold model" , S = V + A, is illustrated along with the 
region consistent with the "signal in noise model". 


4 . 3 Experimental Description 
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Magnitude of Acceleration 
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Figure 4 . 4 Stimulus Diagram and Threshold Predictions 
of Hypothesized Models 
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The experiment best suited to distinguish between 
the "simple threshold" and "signal in noise" models is now 
clear. Stated simply the experiment consists of finding the 
velocity step and acceleration step thresholds for a group 
of subjects and then combining these two stimuli to find 
the subjects’ threshold along the dotted line illustrated 
in Figure 4.4. The predictions of the two models should 
then be compared to the experimental data. 

The experiment was carried out at the NASA Research 
Center at Langely Field, Virginia using a six degree of 
freedom simulator called the real-time dynamic simulator (RDS) . 
All rotations were restricted to the simulator's yaw axis 
which was aligned with the earth' s' vertical axis. The subject 
was seated so that, the axis of rotation was through his head 
and the orientation of his head was varied as necessary to 
achieve either a rolling or yawing stimulus in head axes. 

An analog computer was programmed to rotate the RDS in yaw 
with either a velocity step or an acceleration step whose 
magnitude and onset time could be controlled by the experi- 
menter. Position and velocity feedback from the simulator 
were recorded throughout the experiment and enabled estima- 
tion of actual stimulus magnitude with an accuracy of better 
than + 2/£. Two way voice communication with each subject 
was available in addition to two hand held switches used 
for subject responses. 
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Five men and one woman served as subjects. All attested 
to their good health and normal hearing . Vestibular function 
was tested by means of the sensitized Rhomberg test in which 
subjects must maintc.in their balance for one minute while 
standing toe to heel with their eyes closed. 

For any one experimental session the subject’s head posi- 
tion and the stimulus type (velocity step or acceleration 
step) was held constant. Experimental stimuli were presented 
following the random double staircase method, described by 
Clark and Stewart (Ref. 15 ) except that subjects were in- 
structed to refrain from making responses for cases in which 
no positive sensation of rotation was felt and for which 
the response would be nothing more than a pure guess. In 
such cases the lack of a response was interpreted as if an 
incorrect response had been given. To aid in maintaining 
subject alertness , the subjects were asked to give a graded 
response which would reflect the confidence they had in the 
correctness of that response. If the subject was sure that 
he was rotating to the right (left) he would depress the 
switch in his right (left) hand twice in quick succession. 

If he was moderately confident that he was rotating to the 
right (left) he would depress his right (left) switch only 
once. If the subject was unsure, but had some basis for 
believing his motion was to the right (left ) , he would de- 
press first his right (left) and then his left (right) switch. 
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One experimental session consisted of approximately 
70 stimulus presentations with a short break after the first 
40 were completed, and lasted slightly over an hour. All 
subjects wore light tight goggles and were instructed to 
keep their eyes open during the experiment. Each presenta- 
tion began with the experimenter announcing that an experi- 
mental motion would begin within 15 seconds followed by a 
random stimulus onset delay of between zero and 15 seconds. 
The subject was allowed 13,9 seconds to make a response and 
the first response was the only one accepted. In a very few 
cases a subject reported that he had indicated opposite to 
that which he intended and this was generally accepted by 
the experimenter and corrected for. 

After each presentation the simulator was returned to 
its zero position (in most cases the subject could detect 
this return and thus surmise the direction of his last 
motion) after which followed a period of no motion lasting 
at least 30 seconds. 

The set of all stimulus presentations following the 
point at which the staircases either met or crossed served 
as the basis for estimates of thresholds. 

Figure 4.5 shows a typical strip chart recording for 
a velocity step stimulus in which the subject, while unsure 
of his decision, indicated 2 seconds after onset that he 
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was rotating to the left. Although the tachometer output 
was quite noisy, because of the small rotation rates, it is 
clear that the response of the simulator was quite good. 

The actual rotation rates were calculated based upon the 
information fed back from the position potentiometers. 


4 . 4 Analysis of Experimental Results 

Figure 4.6 shows the thresholds for each of the six 
subjects. Thresholds were calculated as follows: 
Estimated threshold* 3 


TH 


10 



(4.18) 


Upper Confidence Limit = 


N 


£ (lo9 10 (Th)-log 10 s i ) /N 
i=0 


TH X 10 

Lower Confidence Limit * 


= THa 


(4.19) 


N 2 

I (log 10 (TH)-log 10 s. ) /N 
i=0 


TH X 10 


= TH/a 


(4.20) 


where = stimulus magnitude of i th stimulus since stair- 

cases met or crossed 

N+l = number of stimuli magnitudes s Q , s^, 

The value for N differed for each session since it 
dex>ended on how rapidly the staircases crossed but usually 
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Figure 4.6 Threshold Data for Sis* Subjects 
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it varied between 29 and 32. 

Because of Weber's law {Ref. 31 ) all statistics were 
based on the log of stimulus magnitudes. For this reason 
the upper limit on the estimate of thresholds should be per- 
ceptually just as much greater than the estimated threshold 
level as the estimated threshold level is greater than the 
lower limit. 

The number printed next to each threshold estimate 
represents the experiment number for that subject. Since 
combination stimuli for each subject could not be used until 
their thresholds were determined for velocity and accelera- 
tion steps, it was impossible to design the experimental 
presentation in such a way as to rule out order effects. 

To check for learning, each subject's velocity step threshold 
was tested several (2 or 3) times. The first velocity step 
experiment was always either the first or second time in the 
simulator. The remaining velocity step tests were scattered 
among the remaining tests, so that by using the results of 
all six subjects, a model for the effects of learning could 
be formulated. The assumption is made that all subjects had 
approximately the same learning profile and that this profile 
affected all experiments uniformly. 

The model for learning used to correct the experimental 
data is given by: 

L (N) = [K + <1 - (4.21) 

where the threshold resulting from the w session should 



100 


be L(N) times the threshold which would have been found if 
the same experiment had been done during the first session . 
Notice that after an infinite number of sessions (M = «) the 
threshold should. be K times the initial threshold. The 
parameters K and M were varied to minimise the weighted 
variation of the models prediction from the experimental 
data. The optimal values of K and M were found to be 

K* = ,52 and M* ® 3,32 (4,22) 

Figure 4.7 shows a plot of Xog^ 0 (L* (Kf ) ) with the rele- 
vant estimated thresholds and confidence limits shown. Note 
that the confidence limits are symmetrically placed with 
respect to the estimated thresholds when plotted logarithm- 
ically » 

It should be mentioned here that the results of one 
experimental session among the 35 conducted was rejected. 
Specifically, in the second session with subject number four, 
the staircases wandered to extremely high values just before 
the break and when questioned about it the subject admitted 
that he had nearly fallen asleep. After the break, the 
staircases immediately plummeted to reasonable levels, but 
the wide variation in magnitudes (note the confidence limits 
for this threshold) coupled with the subjects' remarks indi- 
cated that it would be reasonable to disregard this session. 

The function (N) can now be used to adjust all 
threshold estimates to reflect the threshold which one would 




expect to find if the subject had been in the simulator for 
the first time* The data could have been adjusted to reflect 
any degree of experience,, but it was felt that most cues 
which contribute to learning in these experiments are non- 
vestibular (tactile, audio, stc.) . Table 4.1 summarizes 
the estimated thresholds and their corrected values . Weighted 
averages for each experiment type wer© calculated and the 
results are as follows (all motions are about an earth vert- 
ical axis ) i 

Average - velocity step threshold for 1.1 . ! = oSSSVsec 

2 

Average acceleration step threshold for Roll = .190°/sec 

2 

Average acceleration step threshold for Yaw = J60°/sec" 

The threshold for yaw was computed based on only five 

of the six subjects since subject number three was not tested » 

If a direct comparison is to be made between the vestibular 

sensitivities to acceleration steps about the yaw and roll 

axes the threshold for the roll axis should be calculated 

based upon the same five subjects . When this is done, we 

2 

obtain a roll acceleration threshold of 0 l58°/sec „ Clearly 
there is no significant difference between this result and 
the . 160°/sec found for the yaw axis= This result supports 
conclusions made by Clark and Stewart (Ref 0 16 ) that there 
are "no statistically significant differences among the 
thresholds about the three major body axes 6 ’* 
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TABLE 4.1 


SUBJECT 

NUMBER 

EXP 

TYPE 

THRESHOLD , 
DEG/SECU-EG/SEC* 4 ) 

CORRECTED FOR 
LEARNING 

0 

FACTOR 

LOG 10 ( FACTOR) 


V 

rm 


1.419 

.15192 


A 

(.151) 

(.173) 

1,445 

.16026 


c 

,193 (. 054) 

. 247 ( . 069) 

1.585 

.20030 

1 

V 

.498 

.697 

1,549 

.18959 


Y 

(.157) 

(.236) 

1.486 

.17186 


V 

.181 

.289 

1.724 

.23665 


V 

.906 

.906 

1.479 

.17029 


A 

(.139) 

(.159) 

1.393 

.14388 


C 

.46, (.109) 

.59 (.139) 

1.531 

.18492 

2 


.30 (.044) 

.42 (.062) 

1.660 

.21973 


a 

.527 

.794 

1.361 

.13372 


H 

(.042) 

(.067) 

1.387 

.14235 


V 

.455 

.760 

1.456 

.16374 


A 

(.299) 

(.299) 

1,596 

.20298 


V 

1.007 

1.151 

1.432 

.15995 

' 

3 

c 

.75 (.237) 

.96 (.303) 

1.367 

.13577 


V 

.780 

1.09 

1.375 

.13843 


A 

(.129) 

(.129) 

1.743 

.24739 


V 

. 941 

1.075 

1.786 

.25203 


V 

.337 

.431 

1.574 

.19692 

4 

. 

c 

.22 (.074) 

.31 (.104) 

1.545 

.18890 


Y 

(.107) 

(.161) 

1.276 

.10605 


V 

.374 

.597 

1.485 

.17166 


V 

.800 

.000 

1.491 

.17352 


A 

(.111) 

(.127) 

1.439 

.15814 


c 

.29 (.070) . 

.37 (.089) 

1.380 

.14003 

5 

Y 

(.131) 

(.183) 

1.446 

.16033 


C 

.63 (.081) 

.95 (.122) 

1.400 

.14624 


V 

.670 

1.069 

1.236 

.09210 


A 

U279) 

(.279) 

3 .596 

.20315 


V 

.957 

3 .094 

1.283 

.10833 


C 

.61 (.009) 

.65 (.114) 

1.675 

.22393 

G 

Y 

(.125) 

(.175) 

2.061 

.31411 


C 

.30 (.095) 

.45 (.143) 

1.429 

.15517 


V 

.500 

.926 

1.337 

.12607 


KEY t 


V Roll Velocity Step 
A Roll Acceleration Stop 


C Combination Roll 
Y Yaw Acceleration step 
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"signal in noise" models of central processing , the data is 
first corrected for learning and then transformed so that 
when plotted on a single stimulus graph (such as shown in 
figure 4*4) all velocity step thresholds fall at and all 

acceleration step thresholds fall at The data from the 

combination stimuli can then be directly compared to the 
model predictions 'illustrated in figure 4*4* 

Figure 4*8 shows the combination stimuli with their con- 
fidence limits plotted- as described above.* The predicted 
thresholds for the two models are shown (an intermediate 
value of K = . 6 is used for the "signal in noise model") 
along with a curve marking the midpoints between the two 
predictions (the decision boundary) * 

Each estimated threshold point plotted can be considered 
the mean of a probability distribution for the true threshold 
value with the confidence limits being the one standard 
deviation points. 

Since this graph shows only relative variations between 
experimental results, any feature of the experimental apparatus 
procedure or statistical technique which causes all results 
to be affected by a given scale factor will not affect the 
validity of any conclusions drawn from the graph* 

If the errors which give rise to these distributions 
are independent (except for an arbitrary scale factor , which 
as noted above, will not affect the conclusions) and gaussian. 
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Figure 4,8 Comparison of Data to Thresholds 
Predicted by Hypotheses 
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then the relative likelihood of the two models being the 
source of these experimental results can be easily calculated. 

One method of choosing between the two models is to ask 
the question: What is the probability that the true thresholds 

are inside the decision boundary? Using the above mentioned 
assumptions, the probability that the true thresholds are 
inside the decision boundary (favoring the "signal in noise 
model") is 99. 8S. 

It is clear from this calculation that if a choice is 
to be made between the proposed models, the "signal in noise 
model" is to be much preferred over the "simple threshold 
model" , 

4 . 5 Conclusions 

The experiments described in this chapter yielded the 
following threshold informations 

1. Roll velocity step threshold « .335°/sec 

2 

2. Roll acceleration step threshold - .19Q°/sec“ 

2 

3. Yaw acceleration step threshold » «160°/sec (4.23) 

and supported the following conclusions: 

1. There is no significant difference in humans between 

sensitivity to rotation in roll and sensitivity to 

rotation in yaw. 

2. The phenomenon of vestibular thresholds can be 

accounted for by a model of central processing of 
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vestibular information consisting only of an optimal 
processing of afferent firing rates in additive noise 
with no necessity for peripheral dead zone nonlinearities . 
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Chapter V 

Stochastic Model for Detection of Hear Threshold Rotational 

Stimuli 

Once the "signal in noise" hypothesis is accepted as an 
adequate model for the basic mechanism of the vestibular thres- 
hold phenomenon * it is possible to hypothesis® a stochastic model 
of vestibular perception which is valid for near threshold stim- 
uli* In the course of creating this model* several significant 
problems arise among which are the followings 

X« Based on the response of the afferent model for the 
semicircular canals developed in Chapter Thr.ee* one would 
predict that the threshold level for a velocity step stim- 
ulus would be between 10 and 15 times greater than the 
threshold level for a step in acceleration • The experi- 
mental results described in Chapter Pour indicate an average 
factor of between 4 and 5 (the ratio of velocity to acceler- 
ation step thresholds for subjects taJ^en separately weres 
2*85* 5«14* 3.94* 4.03* 7.68* and 3 0 64) 0 This discrepancy 
between threshold and suprathreshold perception must be 
predicted by our model of vestibular perception* 

2*. While the peak afferent response for a step in angular 
acceleration occurs approximately 20 seconds after stimulus 
onset and does not drop to half peak value until 45 seconds 
after onset* maximum detection latencies are almost always 
less than 15 seconds (Ref *13 *36 ) . Therefore* some mechanism 
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must be found to explain why average threshold latencies 
do not reflect the time course of the afferent signal to 
noise ratio. 

3. A mathematical model of the detection processes that 
subjects perform during threshold experiments when they 
attempt to choose between right and left moving stimuli 
must be formulated. It is important that a distinction 
be made between the process of detection and the problem 
of estimating stimulus magnitudes which goes on during 
exposures to larger stimuli. The response of this model 
should exhibit the same basic characteristic® observed 
during threshold experiments and be capable of predicting 
response statistics as a function of time for arbitrary 
near threshold stimuli. 

Section 5.1 deals with the first two problems r while sec- 
tions 5,2 and 5.3 deal with the third problem. Section 5.4 sum- 
marizes the predictions which result from this model of 
threshold detection. 

5.1. First Order Processing 

If the acceleration and velocity step threshold stimuli 
found experimentally (.835 degree/second and .19 degree/second ) 
are applied to the input of the dynamic model for the afferent 
response of the semicircular canals, then the responses shown in 
Figure 5.1 will be obtained. It is perfectly clear that no reas- 
onable detector should find it equally difficult to detect the 
presence of these two signals (as it should, since both are thresh 
old stimuli). Furthermore, the maximum detection latency for a 
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threshold acceleration step using these responses should be 20 
seconds or greater, since it takes that long for the instantaneous 
signal to noise ratio to reach its peak. One relatively 
simple hypothesis which has a direct physical interpretation and 
provides an explanation for these discrepancies was found. 

It is well known that in the presence of no input, there is 
a spontaneous firing rate of the vestibular afferent (shown as 
SFR in Figure 3.1) which must be subtracted before processing by 
the higher centers. Goldberg and Fernandez (Ref. 27 ) found a 
spontaneous rate of approximately 90 firings/second in the squir- 
rel monkey, but the exact value is not important to the argument 
which follows. It is reasonable to assume that this resting dis- 
charge is. not perfectly stable but has some low frequency drift 
associated with it, which, if large enough and allowed to remain 
unfiltered, will give rise to spurious sensations of motion and 
will contribute to the masking of low level stimuli. In the pro- 
cess of eliminating the spontaneous component of firing, the eff- 
ect of low frequency variations can also be reduced by estimating 
them and subtracting them from the afferent signal. Figure 5.2 
shows the proposed pathway which estimates both the spontaneous 
discharge (SFR) and hypothesized low frequency variations (c) and 
subtracts them from the afferent signal. If c and the measurement 
noise (n) are members of stationary Gaussiap random ensembles and 
the signal A equals zero, then the minimum-mean-squared-error 
(MMSK) estimate of c is given by the output of a linear filter. 

It is this filter, (shown between E and F in Figure 5.2) which 



Stimulus 


Afferent 

Dynamics 


Afferent 
A + Firing Ra 




SFR t 


(A?R) 


n 


Low Frequency Drift 
in Spontaneous 


Figure 5.2 First 


First Order Processor 


t \ 



Signal 

Available 

for 

Detection 



113 


accounts for the equalization of the signals shown in Figure 5.1 
and which reduces the spurious sensations of motion which would 
occur if the signal c were processed as if it arose from a real 
stimulus. The assumption that the signal at A is small is, of 
course, not always true so some limitation must be placed on the 
filtering so as not to effect the processing of sUprathreshold 
stimuli. Specifically, if the afferent firing rate (AFR) stays 
within ± 4 a (c referring to the effective standard deviation of , 
c + n) of the resting discharge rate then AFR is filtered as if 
A were zero, but if AFR exceeds these bounds, it is presumed 
that a sufficient signal to noise ratio exists to obviate the 
need for eliminating c. A region comprising ± 4a was chosen 
since in the no stimulus case this region would, for all practi- 
cal purposes, contain c + n 100% of the time. This limitation is 
modeled by the saturation nonlinearity which is shown between AFR 
and E . 

The following procedure is used to specify the statistics 
of c and n and the resulting processor (which includes the satur- 
ation nonlinearity and the optimum processor H°(s))r 

1. For the purpose of maintaining simplicity the process 
which produces c is assumed to be such that c is exponentially 
correlated. This can be modelled as a white process through 
a first order shaping filter as follows; 

Cy 

C( S) = ( c ■ s j 1 )W(s) (5.1) 

2 
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where 


♦^(t) « E[W(t)W{t+xJ] = 6(t) (5.2) 

and thus 

c 2 

4 ca (t) = E[c(t)c(t+T)] =2c^ e "' X ^ /C2 (5.3) 

The measurement noise , n, is presumed to be white and there- 
fores 


* nn (T] = E [n (t) n (t+T ) ] = N 2 6 (t) (5.4) 

Since the optimum processor H° (s) is determined as a function 
of the statistics of c(t) and n(t) we must ultimately deter- 
mine only three parameters and No 

2- If the signal available for processing were c-J-n and the 
object were to produce an estimate of c with minimum squared 
error , then the optimal estimator would be given by the fol- 
lowing wiener filters 


who re 


2 


H°(s5 = 


(°l/ Ne 2> l 
n 2 *n/c 2 ) ( £- +1) 


-2(1 + (Zl) 2, x/2 

C 2 v 'H 


(5.5) 


(5,6) 


Since the actual signal available for processing is SFR + c 
+ n and the object is to estimate (SFR +■ c) we must use a 
filter which passes a constant signal with unity gain « Thus 
the expression for H°(s) given in equation 5-5 is acceptable 
if 

(c 1 /Nc 2 ) 

(n 2 +n/c 2 ) 


1 


(5-7) 
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Actually, the optimal filter gain only approaches unity as 
Cf/N •+ Although unity gain must be used in the ultimate 
model, reasonable values for c^ r c 2 , and N are obtained if 
equation 5.7 is replaced by 

(c./Nc ,) 1 2 

i S— ~ .95 (5.8) 

(n 2 + n/c 2 ) 

The consequence of making this change is to produce a system 
in which SFR is estimated and eliminated with no error while 
the optimal filter gain for estimating c(t) is in error by 
only about 5%. 

3. If the response of 

1 ‘ H ° tsl - l£rf (5.9) 

to- the signals shown in Figure 5.1 are designated s^(t) and 
s 2 (t), then ry is chosen such that 
14 14 

/ s 1 (t)|s 1 (t} |dt - / s 9 (t) |s 9 (t) | dt (5.10) 

0 1 . 0 z 4 

This insures that for a given level of measurement noise and 
spontaneous discharge variation a .835 degree/second velocity 
step and a .19 degree/second acceleration step will be 
detected with equal likelihood during a fouteen second pres- 
entation as was found experimentally. The value of r> nec- 
essary to satisfy equation 5.10 was found to be 
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n - 57 Cl +' Cc x /n} 2 ) 1/2 « 1/3 ( 5 . 11 ) 

4. Since the signals shown in figure Sol are threshold 
signals and are thus correctly detested 75$ of the time* the 
net effect of passing them through the filter {1- E° (s) } 
should be to produce signals s^(t) and ^(fc) which when pro- 
cessed by the detector should yield 75% correct responses. 

The net effect of this condition is to adjust the level of 
c(t) and n(t) together so ap to produce the required level 
of signal masking necessary to predict 75% performance for 
the above stimuli. The detector used for this step is 
developed in sections 5.2 and 5.3. 

5. Equations 5.8 and 5.11 along with the procedure described 
in 4 yield three conditions for ©2 and K. If 5.6 is 
substituted into 5»8 we find that 

2 

{— i) 

- ■■ — — — — « .95 ( 5 . 12 ) 

n * 2. . .. . , c l. 2.1/2 

[1 + ($p) J + U + (jy=-) J 

or that 

c 1 /N - 20 (5.13) 

If the procedure described in 4 is now carried out for 
different values of N while picking c, and c~, to satisfy 5.13 

Jm 

and 5.11 then 75% performance is reached for one afferent 
channel when 

c^ = 4.46 

c 2 = 60 (5.14) 

N - . 223 
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Since n = 1/3, the optimal filter is given by 

H °< s > = 3HTX (5 * 15) 

The linear region for the saturation nonlinearity is bounded 
bySFK±4a where a is the effective standard deviation of the 
signal c(t) + n(t). Since 

c 2 

4 a = 4(^i- + N 2 ) 1//2 = 2 ips (5 - 16) 

2 

The region over which the optimal processor is allowed to 
operate freely is given by fcFR-4a, SFR + 4a) = (88, 92) ips. 
Figure 5.3 summarizes the slowly varying spontaneous dis- 
charge and the first order processing. It is interesting to note 
that while this system was designed to meet several key criteria, 
it also meets several conditions which although, not specifically 
imposed, certainly would be expected from acceptable models of 
vestibular processing. 

For example: 

1. It might have turned out that the reasonable region 
for signal processing (90±4a) would not have been 
sufficient to include the maximum variation in firing 
due to threshold stimuli. Inspection of Figure 5.1 
shows that this is not the case since the peak vari- 
ation in firing due to a stimulus of 0.19 degree/ 

2 

second is 0.78 firings/second which is well within 
the specified region of ±2.0 firings/second. 
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Figure 5.3 Model Of Information Available to Detector 
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2. On the other hand, the necessary variation in steady 
state firing is significantly smaller than that found 
in the most regular cells recorded by Goldberg and 
Fernandez (Ref. 27 ). This would lead one to the con- 
clusion that the higher centers must be capable of 
combining the information from several afferent cells 
to achieve the threshold performance found experi- 
mentally. A discussion of this phenomenon is given 
in section 3.1. 

The existence of this three second filter would ex- 
plain why vestibular threshold experiments need to 
last only ten to fifteen seconds. As noted previous- 
ly , if such a mechanism did not exist, the maximum 
signal level from a step in angular acceleration would 
occur 2d seconds after the stimulus onset and not 
fall to half peak amplitude until 45 seconds after 
onset (Figure 5.1). Figure 5.4 shows the signals 

available to the detector for a .835 degree/second 

2 

velocity step and a .19 degree/second acceleration 
step after the first order processing. Not only is 
it believable that these two signals are equally easy 
to detect (as opposed to those illustrated in figure 
5.1) but it is also easy to understand why latencies 
which significantly exceed fifteen seconds are rarely 


encountered. 
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Figure 5 A Signals Available for Detection After 
First Order Processing for Threshold 
Steps in Angular Velocity and Acceleration 
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5.2 Detection Using Information From The Suprathreshold Optimal 
Estimator 

Since the first order processor insures that a threshold . 

velocity step (.835 degree/second) and a threshold acceleration 

2 

step (.19 degree/second ) will be equally easy to detect, and 
since a model consistent with the "signal in noise" hypothesis 
insures that a combination stimuli given by 

w (t) = K ( ,835 + . 19t) (5.17) 

will only yield equivalent performance if 

1/2 < K < 1//? (5.18) 

we can conclude that any reasonable detector which utilizes the 
available signals (G in Figure 5. 3) should meet all of the essen- 
tial requirements. In Chapter Three a Kalman Filter was proposed 
as a model for the higher centers' perception of suprathreshold 
stimuli. Since such a processor is presumed to be operating any- 
way, it would make sense that it could also be used as part of the 
detector for near threshold signals. 

If one looks at the output of the Kalman Filter associated 
with the semicircular canals at any given time it is seen that it 
consists of an estimate of the rotation rate and an error variance 
for the estimate. Since all processes are assumed to be Gaussian, 
a probability distribution for the actual stimulus magnitude con- 
ditioned on previous measurements of afferent firing rate is known. 
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At that instant of time* the subject could report on the probabi- 
lity that his instantaneous rotation is to the right or the left 
as illustrated in Figure 5*5- All that is needed is & Hay to pro- 
duce such an estimate under the known conditions of a threshold 

experiment and to. give a model of the subject a s criteria for 
issuing a report* 

In most vestibular theshold experiments a subject assumes 
that any motion given him will either be always .to th© right or 
always to the left and not a combination of the two* -since he is 
expected to only give one report* if the errors in the estimate 
of rotation were independent from one interval to the next* it 
would be an easy matter to combine these estimates to produce the 
probability that the stimulus was to the right (or left) , but 
since these errors exhibit significant correlation from one inter- 
val to the next* combining these estimates properly is a diffi- 
cult computational task* In fact to combine N. such ©©timafee® 
would involve the computation of two N dimensional integrals of 
the joint probability density for the stimulus history conditioned 
on the measurement history. While such a computation could be car- 
ried out by a computer* it seems unlikely that the brain would be 
doing anything similar* or that it would be necessary to use such 
a complicated detector for our model. To avoid this computational 
problem* the Kalman filter is applied sequentially to each meas- 
urement of afferent firing as if it were the only measurement a- 
variable. In this manner we obtain a probability density for w(t) 
(similar to the one shown in Figure 5.5) conditioned on only one 
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Figure 5.5 Decomposition of Probability Density Function 
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measurement at a time* Although this does not reduce the cor- 
relation of estimate errors from one interval zo the next to 
zero, it does eliminate the major contribution to this correla- 
tion (the fact that in the usual Kalman filter, the estimate 
error at one time is propagated and becomes part of the error 
of subsequent estimates ) e With the error correlation signifi- 
cantly reduced the assumption of independence is reinstated and 
the estimates can easily be combined. If . (t.) is the 

probability at time t^ that the motion is to the right given the 
measurement of G at then the combined probability that the 

motion is and has been to the right at time is given by: 


VV = 

and to the left is 


A 

A+B 


(5.19) 


P (t ) SS ■■ 

l'tt a+b 


“ P R (t N } 


(5.20) 


where 


N 

A = i=l PR /G(t i ) (t i ! 


(5.21) 


B= J 1 P L/G(t i ) (t i ) “ < 5 ' 22) 

and P , , (t- ) is obtained by calculating the appropriate area 

R/ G l C ^ ) X 

under as *H ustratec * in Figure 5.5. 

Alternatively (see Ref. 7 ) the log of the ratio of P R (t n ) 
to P L (t n ) can be updated sequentially as follows: 
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In C 


p^rtr 1 


p n-l 


+ 


P R/G (t ) 

ln( 3 — ) 

L/G Ct ) 


(5.23) 


= p n _ 1 + ln(A(G(t n ))) 

P G(t )/R 

where A(G(t )) is the likelihood ratio defined as = . 

n P G(t n )/L 


The above scheme leads to an accumulated probability that the 
rotation .was to the left or to the right, but still does not tell 
us what criteria a subject will use for making a report. 

In the experiments conducted at Langley Research Center and 
reported in Chapter Four the subjects were requested to report 
their direction of rotation only if they had sensations of motion 
significant enough so that their report would not be merely a 
guess. For this reason the minimum criteria for issuing a report 
consisted of P R (t N ) either exceeding .75 or diminishing below 
.25 (implying that P L (t N ) > .75),. Another consideration which 
should be used in choosing a reporting criteria is that subjects 
are aware that they have at least 13 - 15 seconds to issue their 
report. Thus, before 13 seconds have elapsed the criteria for 
issuing a report is stricter than 75% due to the fact that a 
reasonably motivated subject would not report with only 75% con- 
fidence early in the experiment since he would know that his period 
for reporting had significant time remaining and he might as well 
obtain more data before committing himself.- As a result of these 
considerations the reporting criteria shown in figure 5.6 was used. 
It should be made clear that any reasonable decision boundary 
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could be used since the basic integrity of the model is not depen 

dent on this particular profile. The level of measurement noise 

and base rate variation was adjusted, in the manner described in 

2 

the previous section, until a .19 degree/second acceleration 
step was correctly identified approximately 75% of the time. 

The necessary values for and N to obtain 75% performance are 
given in 5.14. As expected, equivalent performance for a 
.835 degree/second velocity step and a combination stimuli 
(W (t) = .67 (.835 + . 19t) degree/second) was reached for the 

w 

same noise level. A typical correct response simulation and a 
typical incorrect response trial are shown in Figure 5.7, 

Due to the numerical complexity of this detector which 
requires the simulation of a Kalman filter, large scale Monte 
Carlo simulations are impractical. Thus, while extensive 
latency histograms as a function of stimulus magnitude and type 
are not available it is possible from the number of simulations 
made to indicate the general nature of the response latencies. 
First it is clear that as stimulus magnitudes are increased, 
latencies uniformly decrease independent of the stimulus type. 
Secondly, the shortest latencies are associated with velocity 
step stimuli, the longest with acceleration step stimuli and 
intermediate latencies are found for combination stimuli. All 
of these predictions are consistent with the available experi- 
mental data for rotational thresholds. 



5. 3 Simplified Detector Model 


The only major drawback of the detection model described in 
Section 5.2 is its numerical complexity and the resultant ex- 
pense required to make a sufficient number of simulations to get 
meaningful statistics. For this reason an attempt was made to 
produce a model which gave equivalent performance* but which had 
sufficient simplicity to make large Monte Carlo, runs realistic. 

The best way to reduce the numerical complexity of the 
detector is to eliminate the Kalman filter and to fry and use 
the available signal in a simpler manner to generate a decision 
parameter „ Since any such simplified detector is likely to be 
less efficient in detecting the stimulus than the. suboptimal 
Kalman detector, we can expect one of two consequences * Either 
1) if the same noise level on the afferent signal is main- 
tained then xnc>re afferent channels will have to be 
processed if the name performance 1© to be achieved, 
or 2} if only one afferent channel is to be used then similar 
performance will require a reduction in the afferent 
noise level. 

These consequences are not significant though, since our object • 
is to obtain equivalent performance with a minimum of computation 
al effort, regardless of the details of the model (e.g» noise 
level, number of afferent channels, etc.) 

Typically, a detector consists of a decision parameter 
(which is a function of the available signal) and a set of 
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reporting criteria which the decision parameter must meet before 
a report is issued. One possible decision parameter could be 
formed by integrating the available signal and issuing a report 
if the parameter exceeded certain bounds. There are two problems 
with this solution. First,- the result of integrating a zero mean 
random process is to produce an output similar to a random walk 
whose variance increases without bound as time increases. Such a 
model would predict a large number of reports in the absence of 
any stimulus input. While such reports are observed experimentally 
they are relatively rare. Secondly the parameter formed by simple 
integration does not take sufficient advantage of the initial 
prominent response to a velocity step stimulus. In fact, its 
expected value asymptotically approaches zero as t goes to infin- 
ity. Inspection of Figure 5.4 reveals that for a velocity step 
stimulus the available signal is of the correct sign and large 
for a short time and of the wrong sign and small for a much 
longer time. The optimum Bayes detector for the case in which 
one wishes to distinguish between two uncorrelated signals with 
different variances consists of a square law device and an integ- 
rator (Ref 69, 71 ). Such a detector tends to amplify the 
importance of large deviations more than small ones. This charac- 
teristic is useful in the case of a velocity step since it weighs 
the initial prominent signal much more than the small undershoot. 
One drawback of this procedure is that by squaring the available 
signal the sign of the signal is lost and thus it would be impos- 
sible to distinguish between right and left moving stimuli. 



The decision parameter must therefore meet four key criteria g 

1, It must be able to distinguish right moving from left 
moving stimuli* 

2 a it must have a bounded variance in response to a 
no-stimulus trial to reduce the probability of a 
false alarm to a level consistent with experimental 
data* 

3, The processor which gives rise to the decision para- 
meter must have a form consistent with that analysed 
in Section 4c 2 and therefore be consistent with the 
M signal in noise model" predictions illustrated in 
Figure 4,4, 

4 a It should incorporate some technique for taking 
advantage of the larger magnitude of the initial 
response as compared to any subsequent undershoot, 

A detector which satisfies these four conditions and which 
is numerically simple is shown in Figure 5*8. The first condition 
is satisfied since the product G | G(t^) | which drives the 
processor which generates H retains the same sign as G(t<), The 
variance of H is bounded in response to random processes with 
bounded variance since H is the output of a discrete time first 
order filter, A filter time constant of four seconds (e ~°78) 
was chosen to reduce the false alarm rate to less than 10%. Thus 
decay mechanism also serves as a memory limiter since it progres- 
sively discounts the importance of information as that information 
grows older. The third condition is fulfilled since the linear 



H(t k )=.78H(t k _ 1 ) 

+G(t k > |G(t k ) 
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operator which processes G(t} jG(t) | satisfies the requirements 
of a "signal in noise" processor as given in Section 4o2. Fin- 
ally , since G(t) 5 s contribution to a change in H is proportional, 
to the square of the magnitude of G(t), H responds with increasing 
sensitivity as | G ( t ) | increases* 

The minimum computational effort will be reached if only one 
channel is simulated. For this case (N=l) the measurement noise 
and resting discharge variation must be reduced fey a factor of 
^ /I so that for this simplified detector 


= 2,58 
N = o 129 


(5 o 24) 


and the equivalent reduction capability of the central processor 
would be given by (see Section 3*1.2) 

Reduction Capability = CTT2T 5 2 ° 1560 (5.25) 

which is three times greater than was necessary for the suboptimal 
detector using the Kalman filter. The response delay has not 
been modelled but it is expected that the delay is longer than 
typical reaction times since it arises from the conscious weigh- 
ing of barely noticeable perceptions. In fact, the delay may be ! 
a function of how much the detection parameter, H< exceeds the 
decision boundary since a large H would imply an easier decision 
and thus less deliberation. In any case a Rayleigh or Maxwell 
distribution for the response delay with a mean delay of one or 
two seconds should yield reasonable predictions concerning re- 
soonse latencies, , Finally the decision boundary (D) was adjusted 
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2 

to insure that an acceleration step of 0.19 degree/second 
would be detected 75% of the time. 

The model shown in Figure 5.8 was simulated using the 
parameter specifications outlined above. The computational 
simplicity is revealed by the fact that over two thousand, 
twenty second simulations could be run with approximately one 
minute of central processor time. The results of simulations 
using various stimuli of differing magnitudes is summarized in 
the next section. 



5 « 4 Summary of Model Predictions for Detection Probabilities 
and Latencies 


Velocity step, acceleration step and combination stimuli 
were tested to determine the model 8 s predictions concerning the 
probability of detection during a fourteen second trial and the 
distribution of detection latencies » Figure 5.9 summarizes the' 
performance predicted by the model as determined by Monte Carlo 
simulation. Each point represents the results of one hundred 
simulated trials « Stimulus magnitudes are referred to the 
experimentally determined threshold levels as follows? 

a (.8.15) deg/sec velocity step (5.26) 

THa N = ^ a N ( • 19t) deg/sec acceleration step 

I (K ( . 835+ « 19t) deg/sec . combination 

where a = 10 = .1259 

K » .67 (determined from Figure 4.12) 
and N = -2, -1,0, 1,2 

While it is impossible to determine the precise threshold 
predictions in this manner it is clear that the predicted thres- 
holds (755. correct response) are between TH /o and THc for all 
throe stimuli. It is unnecessary to determine the exact thres- 
hold predictions with any more precision than this since much 
larger variations are seen experimentally among subjects. 

Figures 5.10, 5.11 and 5.12 illustrate typical time histories 
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of H ( t) for each of the three stimuli* In each case one correct 
response and either a false response or no response trial is 
shown. 

Figures 5.13, 5.14 and 5.15 summarize the model's pre- 

dictions for response latencies at threshold (no reaction delays 
are included in these figures). Short latencies for velocity 
steps, relatively long latencies for acceleration steps and 
intermediate latencies for combination stimuli are typical of 
the results. It is clear from the model that in all cases 
average latencies will decrease as stimulus magnitudes increase. 
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Figure 5„14 Latency Histogram for Threshold Acceleration Step 
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CHAPTER VI 

STOCHASTIC MODEL FOR DETECTION OF NEAR THRESHOLD CHANGES IN 

SPECIFIC FORCE 

The model for threshold detection of rotational stimuli 
was found to consist of three distinct components. The first 
is a linear model of the afferent response of the semicircular 
canals. The purpose of the second component, referred to as 
the first order processor, is to eliminate the spontaneous 
firing rate and any associated low amplitude, low frequency- 
variations from the afferent signal. The third and final 

i.f. 

section models the detection and decision processes which 
utilize, the signal available at the output of the first order 
processor. Detection of changes in specific force can be 
modelled using the same fundamental approach except that the 
relevant sensory signals are now those associated with the 
otolith organs. 

A linear model relating changes in specific force to 
changes in the afferent firing rates from the utricular and 
saccular organs was developed in section 3.2 and will be em- 
ployed here without modification. The requirement for a first 
order processor, similar to that found necessary for the pro- 
cessing of semicircular canal information, will be discussed 
in 6.1. While a suboptimal detector employing the appropriate 
Kalman filter could be developed for the otoliths in a manner' 
similar to that done for the semicircular canals (Section 5.2) 
it was decided that the benefits of this approach were far 



outweighed by the associated computational complexity „ Instead* 
a successful attempt was made to implement a detector similar 
to the "simplified detector" developed in Section 5,3* The 
specifications of this detector are outlined in Section 6.2 
and Section 6.3 summarizes the resulting model predictions. 
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6 . l Necessity for First Order Processing 

The problem of detecting linear acceleration requires 
the higher centers to distinguish changes in the firing rate 
of first order otolith afferents. If, in the absence of 
acceleration or tilt, the discharge level was truly constant . 
and if this level was precisely known by the higher centers 
than it could simply be subtracted and the resultant signal 
processed. On the other hand, if the spontaneous discharge 
level is not stable, but instead varies slowly within some 
known bounded region then some mechanism must be postulated 
to distinguish between small changes in the afferent firing 
caused by stimulation of the otolith organs and those caused 
by normal variations in the spontaneous discharge level. 

There -are several approaches which the system could take to 
solve this problem. One approach would be to ignore the fact 
that some low frequency low amplitude changes in the afferent 
firing rate are not stimulus related and merely subtract the 
average spontaneous discharge rate from- the afferent firing 
rate and therefore process the resulting signal as if it arose 
from a true input stimulus (see Figure 6.1). The problem with 
this approach is that one would continually be alerted to 
accelerations and tilts which did not actually take place. Such 
a high false alarm rate is unacceptable and as in the case of 
"the boy who cried wolf" some mechanism will arise to insure 
that such continual false alarms will be ignored. If these 
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low frequency variations occur, say 99% of the time within 
some bounded region of the average spontaneous firing rate 
then a second approach (see Figure 6.2) would be to ignore 
completely any change in the faring rate which takes place 
within this bounded region. This approach would solve the 
problem of a high false alarm rate since the signal avail- 
able to the detector would only rarely differ from zero when 
there was no corresponding stimulus present. The problem 
with this approach is that too much information is thrown 
away unnecessarily. While it is true that smali changes in 
the afferent firing rate which have the same spectral compo- 
sition as those arising at c (Figure 6.2) but which are due 
to true stimulus inputs cannot £e distinguished from those 
arising at c , it is not true that signals with a different 
spectral composition could not be distinguished. The simp- 
lest model for a spectral decomposition limited to variations 
in the firing rate near the known average rate of spontaneous 
discharge is shown in Figure 6.3 0 The model shown in Figure 
6.3 is actually a compromise between the model shown in Figure 
6.1 and that shown in Figure 6.2. In the first model, the 
presumption is that all variations in the firing rate regard- 
less of their frequency content will be assumed to arise from 
stimulus inputs to the sensory organs. In the second model 
the presumption is that all variations within a bounded region 
regardless of their frequency content will be presumed to 
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arise from the process c and mupt therefore be eliminated. In 
the third and final model the presumption is made that for 
small variations in the firing rate the low frequency portion , 
of the signal variations arise from the process c and the high 
frequency portion of the variations arise from a stimulus 
which must be detected. To this end the first order processor 
passes high frequency variations while eliminating low fre- 
quency variations. 

It is interesting to note that this first order processor, 
while essentially eliminating the direct effects of any low 
frequency variations in the spontaneous discharge; preserves 
the systems ability to detect the most important class of stim- 
ulus disturbances? specifically abrupt changes in specific 
force which will lead most rapidly to a significant change in 
one ' s position . 

The above analysis of first order processors is equally 
applicable to the modelling of the perception of near threshold 
rotational stimuli based upon semicircular canal afferents. For 
the case of rotational stimuli the first orcler processor was 
postulated to account for the unexpected ratio of the velocity 
step threshold to the acceleration step threshold (Section 5.1). 
No such experimental discrepancy exists for the otolith organ. 
The only significant implication of the first order processor 
for threshold detection of steps in specific force would be a 
slight shortening of the predicted latencies. Figure 6.4 
shows the signal available for detection in the case of a 
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igure 6.4 Signal Available for Detection for Various 
Values of t 
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threshold step in specific force (.005 g) for various values 
of the filter time constant, x. Figure 6.5 shows the result 
of passing each of these signals though the detector used for 
the semicircular canals (Figure 5.8, N=l) and normalizing 
their peak responses. Since the average detection latency will 
occur about the time that H(t) peaks it can be concluded that 
as t decreases the average predicted latency will tend to de- 
crease. It is difficult to choose a value for x based on 
these curves since they are only approximate indicators of 
predicted latencies and since the predicted latencies do not 
change radically with changes in x. In fact, given a small 
leeway in choosing the statistical distribution of the response 
delay T, any of the values for x shown would be in reasonable 
agreement with the experimentally determined latencies (Ref. 51 ). 
If the mean response delay were presumed to be between one and 
two seconds for near threshold stimuli then any value of x 
greater than or equal to three seconds would be reasonable. In 
section 6.3 the model predictions will be given for t= 3 seconds 
and t= 30 seconds to show the effegt of x on the predicted lat- • 
encies. 

Specifications for the low frequency variations in the 
spontaneous firing rate can be determined in the same way as 
was done for the semicircular canals in Section 5,1. The equa- 
tions which must be satisfied by N, c^ and c 2 (see Figure 6.3 
and Section 5.1) ares 
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Figure 6.5 
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Normalized Detection Parameter for Various 
Values of t 
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or 


N 


« 20 


( 6 . 2 ) 


and 


[1 + (“H 


c,s2 1/2 
' ] 


1/t 


(6.3) 


N was determined in Section 3.2 to be 
N - .147 

Thus from Equation 6.2 
c x = 2.94 

and from Equation 6.3 


c 2 * 20t 


(6.4) 


(6.5) 


( 6 . 6 ) 


The only parameter which remains unspecified in Figure 
6.3 is the region over which the first order processor oper- 
ates linearly. The linear region for the saturation nonlinear- 
ity shown in Figure 6.3 is bounded by SFR±4c where o is the 
effective standard deviation of the nonstimulus related signal 
variations. a is given by 


a 


2 1 
N + 2^ 


1/2 


(6.7) 
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Sc 2 Simplified Detector Model 

% 

The simplified detector which was used to model the de- 
tection of near threshold rotational stimuli (see Figure 5 0 8) 
is also capable of predicting threshold performance for detec~ 
tion of steps in lateral acceleration „ Sine© the minimum 
computational effort is required when only on© channel is 
simulated, N is set equal to 1* The motivation for the struc- 
ture of this detector is outlined in detail in Section 5.3 
and is not repeated here, 

“Ts . 

The response delay e 9 which represents the sum 
total of all delays which are not related to the time history 
of the signal available for processing* is not modelled in 
detail here, except to remark that T should be. viewed as a 
positive random variable whose minimum value should be set 
equal to three or four hundred milliseconds (which represents 
the time normally associated with supr&threshold reaction 
times) and whose mean value should be set to one or two sec~ 
onds. The value of D p which represents the minimum deviation 
of the decision parameter H(t) which will provoke a response, 
should not differ significantly from the value- found for the 
semicircular canals. The exact value of D can be chosen in 
two ways. Either: 

lo D can be adjusted to yield a threshold prediction 
which is the same as the experimentally determined 
value and then D compared to the value found for the 
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semicircular canals to check, if it is reasonable, 
or 2. D can be set to the valiie found for the semicircular 
canals, the model 9 s theshold prediction determined 
and then this threshold prediction compared to the 
experimentally determined thresholds* 

For the case t=* 3 seconds, which is the same time constant used 
in the first order processor for the semicircular canals, there 
is no need to apply both procedures since for D = 0.16 (the 
value used for the semicircular canals) the resulting threshold 
(75% correct detection) is .005 g for the utricles and .010 g 
for the saccule. These values are in complete agreement with 
those of Meiry (Ref. 51 ) . For x«30 seconds, a value of 
D = .32 is necessary if the above threshold values ate desired. 
This is only a factor of two greater than that used for the 
semicircular canals. If D is set equal to .16 then for t «30 
seconds the predicted utricular threshold would be approximate- 
ly .0035 g and the predicted saccular threshold would be .007 g. 

■ i, . 
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6. 3 Summary of Model Prediction for Detection Probabilities 

and Latencies 

Acceleration steps of five different magnitudes were 
tested to determine the model’s predictions for the probabi- 
lity of detection during simulated fourteen second trials. If 
"TH" is used to designate an acceleration step of oO05 g then 
the utricular model was tested with steps of magnitude ‘£H/o 2 , 
TH/o , Til, THcr and THo 2 where c** 10* \ Figure 6*6 summarises 
the results of these Monte Carlo simulations. For the simu- 
lations in which T“3 seconds the decision level D was set to 
.16 and for t= 30 seconds D was set to .32, The predictions 
are essentially equivalent with 75% performance occurring for 
a stimulus magnitude of .005 g. Each data point shown in 
Figure 6.6 is the result of 210 simulated trials-. Using the 
same value for D in the saccular model will result in essen- 
tially the same curves shown in FigurO S„S esceopt that TH 
would have to be set equal to .010 g sine© the sensitivity of 
the saccular organ is only half' that of the utricular organ 
(see Section 3.2). Thus the threshold prediction for the 
saccules would be. 010 g. 

Figure 6.7 shows the time history of the detection para- 
meter H ( t ) for a correct response trial and an incorrect re- 
sponse trial for the model with t= 3 seconds. Figure 6.8 shows 
the same simulations of the model for t=3Q seconds.. In both 
figures the stimulus magntiude is .005 g. 
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Stimulus Magnitudes (a-lO’ 1 ) 


Figure 6.6 Model Predictions of Performance Variations 
as a Function of Stimulus Magnitudes 


Detection Parameter Detection Parairt 


154 



Figure 6.7 Threshold Acceleration Step Simulations 
( t =3 seconds) 



Figure 6.8 Threshold Acceleration Step Simulations 
(t =30 seconds) 
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Figures 6.9 (t=3 seconds) and 6.10 (x=30 seconds) sum- 
marize the predicted distributions of the response latencies 
due to the time history of the signals available for detection. 
These distributions do not include the response delays modelled 
by e Ts . Inspection of Figures 6> 9 and 6.10 reveal that for a 
tenfold increase in t the predicted increase £n the average 
response latency would only be, approximately one second. Since 
this dif ference . latency^can. j yet easily be incorporated 

in the term e the value qhosen for^T %Ls not very critical. 

One of the most interesting implications of these results 
is that the detection processes which make use of the signals 
available at the afferent level for the semicircular canals 
and the otoliths can be modelled as if they were identical. In 
both cases the afferent signals were corrupted with white 
measurement noise and low frequency variations in the sponta- 
neous discharge. In both cases a Weiner filter was used to 
estimate these low frequency variations which .were then sub- 
tracted from the afferent firing level to produce a signal 
which presumed to consist of a white measurement noise process 
and possibly a signal related to a stimulus acting on the 
sensor. Finally, identical detectors with the same detection ; 
parameters and the same criteria for issuing a report were 
used to predict near threshold detection performance consistent 
with experimental data. 
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LATENCY (SECONDS) . 

Figure 6.9 Latency Histogram for Threshold 

Step in Acceleration (T®3 seconds) 
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Figure 6.10 Latency Histogram for Threshold 

Step in Acceleration (t =30 seconds) 
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CHAPTER VII 

PERCEPTION OF STATIC ORIENTATION IN A 
CONSTANT SPECIFIC FORCE ENVIRONMENT 

In previous chapters the perception of rotation or ac- 
celeration was always presumed to originate from the simple 
stimulation of a synergistic pair of sensors {e.g* the right 
posterior canal and the left superior canal) and no considera- 
tion was given to the problem of integrating the sensations 
arising from more than one synergistic pair (e*g. the hori- 
zontal canals and the right posterior and left superior can- 

v 

als) . For the case of the semicircular canal system the 
simplest solution to this problem would be to postulate that 
since the planes of the three canals' are roughly orthogonal 
that the sensations which arise from each synergistic pair can 
be considered the components of a vector which represents the 
net: sensation of rotation. While it is true that the sensitive 
axes of the canals do not coincide with the conventional roll, 
pitch and. yaw axes to which we tend to relate any. sense of 
rotation, the sensation of rotation relative to the conventional 
body axes can be derived using a simple transformation of 
coordinates. 

A similar approach to the problem of integrating the 
sensations from the utricles and the saccules would yield 
a model which was in agreement with the gross nature of the 
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steady state perception of tilt, If the sensitive axes of the 
otolith organs formed a mutually orthogonal set and the 
steady state sensitivities of these organs and their orienta- 
tions with respect to the head were accurately assessed , then 
an unbiased estimate of the head‘s steady state orientation 
with respect to a constant specific force environment could 
be made. (Specific force or gravito-inerfcial reaction force is 
defined as the difference between the gravitational force 
vector and the vector representing the translational accelera- 
tion with respect to inertial space.} Unfortunately experimen- 
tal data suggests that the perception of tilt is biased and 
that this bias is a function of the intensity of the specific 
force field (lg,2g, etc.) and the orientation of the head 
relative to the local direction of specific force. 

The purpose of this chapter is to discuss the nature of 
these biases and to suggest a simple model for the integration 
of otolith information which will yield predictions consistent 
with the available experimental data. 
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7.1. Perceptual Illusions of Static Orientation 

The first systematic investigation of man's perception 
of the vertical was carried out by Mach [ 45 , 46, 47] 
beginning in 1873. In Mach's experiment the subject was 
rotated in a centrifuge at a constant rate while facing in 
the direction of rotation. The net specific force was such 
that the subject invariably felt tilted away from the center 
of rotation. When a plumb bob was suspended inside the 
centrifuge Mach reported that his sensation of the vertical 
was generally aligned with the direction indicated by the 
plumb bob (i.e. aligned with the local direction of specific 
force) . 

In more recent years numerous psychophysical experiments 
have been conducted, e.g.,[ 11, 12, 29, 30, 52, 56, 57, 

63, 64, 72, 74, 75 ] in which subjects were tilted with 

respect to the local direction of specific force 1 and then 
requested to indicate their perceived vertical or horizontal. 

In addition Cohen [ 17 3 and SChbne [63 3 investigated the 
effect of specific force magnitude and pitch orientation on the 
perceived angle of pitch. 

The mechanism used to generate sustained tilt angles and 
the apparatus employed to permit the subject to indicate his 
orientation are crucial to the interpretation of the experimen- 
tal results. For experimental results to be useful in modelling 
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the static sensation of orientation as a function only of 
the output of the otolith organs several key criteria must 
be met . 

1. There must be no visual information available to 
the subject which in any way affects his percep- 
tion of orientation. The ideal way to meet this 
condition would be to have the subject's eyes closed 
but unfortunately this eliminates the most reliable 
methods for the subject to indicate his orienta- 
tion. Therefore the experimenter must strive to 
design an indicator which has the least influence 

on the subject’s perception of orientation. One 
common solution to this problem,, while probably 
not optimal , is to display a luminous line which 
the subject can manually rotate by means of a knob 
to indicate his perceived vertical in an otherwise 
darkened enclosure. 

2. Data concerning the subject’s perception of orienta- 
tion must be collected after the experimental 
conditions have stabilized and have remained constant 
for a period sufficient to abolish all response 

from the semicircular .canals. While one minute or 
more would be preferred a period of no less than 
thirty seconds should be adequate. There have been 
experiments reported in which a subject seated in a 
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centrifuge is slowly accelerated while instructed 
to continuously report his perceived orientation. 
While the results of three experiments are useful 
indicators of the phenomenon we are investigating, 
they should be analyzed carefully with regards to 
the possible contaminating effects of both canal 
involvement and the dynamic portion of the otolith 
response; especially at times just following abrupt 
changes in acceleration. 

3. Since the effect of proprioceptive information on 
perceived orientation is not the object of study 
here, it is important that this source of information 
be eliminated insofar as possible through suitable 
head and body supports. 

and 

4. Instructions to the subject must clearly indicate 
that his task is to indicate the orientation of the 
earth vertical and not his perception of the direc- 
tion of net force. 

With these criteria met and a clear record of both the 
specific force magnitude and the orientation of the head 
with respect to the local direction of specific force a 
reasonably consistent set of experimental data emerges. 

In a certain class of experiments in which a subject is 
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tilted laterally with respect to the local direction of 
specific force the subject perceives that his tilt angle is 
less than the objectively measured angle* This effect 
is named the Aubert effect (or A effect [3 ])<> The Aubert 
effect manifests itself in several ways= If for example 
a subject is tilted to the right in a 1 g environment while 
viewing a stationary vertical line he will s@nse that the 
lino has rotated counterclockwise and is now tilted to the 
left. Alternately if the subject is able to control the 
orientation of this line he will tend to rotate it in the 
same direction as his own tilt in order to have it appear 
vertical to him. In a separate class of experiments a subject 
will experience a sensation opposite to the one described 
above. An objectively vertical line will appear to tilt in 
the same direction as the subject’s tilt and he will thus 
perceive his tilt to be greater than his actual tilt* This 
effect is known as the Muller qffect (or E effect) and was 
first described by Muller in 1916 [53 ]. Although the E 
effect is normally associated with either a transient sen- 
sation or with experiments in which the magnitude of the 
specific force vector is greater than unity, some investigat- 
ors [ 75 ] have reported a miX«| E effect for small roll 
anqles in a 1 g environment. While the terms "Aubert effect" 
and "Muller effect" may be useful in discussing the outcomes 
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of certain experiments they do not by themselves serve to 
clarify the mechanism which underlies the general phenomenon 
associated with the perception of orientation. 

For reasons which will soon become clear all of the 
experiments with which we are concerned will be divided 
into three basic categories. Figure 7.1 illustrates the 
three categories. -Consider a coordinate system whose y 
axis is aligned with the head's pitch axis (pitch forward 
is positive) , whose z axis is perpendicular to the 
average plane of the utricular macula (i.e., approximately 
0= 25—30 degrees pitched back relative to the vertical when 
the head is in a normal erect position [ 19 ]) and whose x axis 
is defined by a right handed cross product 

ix = iy x iz 

If the steady state -specific force vector for a given 
experimental trial is referred to these axes as follows 

SF - SF i + SF i + SF i (7.2) 

— x-x y-y z — 2 

Where SF , SF and SF are expressed in g's 
x y z ^ 

we can categorize that experimental trial according 
to the following rules: 

If SF > - cos ( 0) then designate as category A 
z 

If SF = - cos ( 0) then designate as category N 

ri • 

If SF < - cos ( 0) then designate as category E. 


(7.3) 
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This categorization will be better understood if it is 
noted that for an experiment with the head erect in 1 g the 
stimulus is of type N. Adding any accelerations in the xy 
plane of the coordinate system just described will not alter 
the stimulus type. The crucial test is whether or not the 
component of specific force in the z direction is more 
negative (type E) or less negative (type A) than in the case 
in which the head is erect in 1 g. 

In reviewing the results of several experiments, each 
trial will be assigned to the appropriate category and the 
direction of the mean bias in perceived orientation relative 
to the net specific force noted. As a preliminary example, 
the experiments by Mach described at the beginning of this 
section are of type N since the subject's head was erect in 
1 g and an additional component of specific force was added 
in the negative y direction. Since this change in specific 
force did not alter the z component of specific force the 

r 

category remains -the same as "head erect 1 g " — namely cate- 
gory N. Mach's experimental result was that the vertical was 
aligned with the net specific force without any net bias. 

More recent centrifuge experiments have confirmed Mach's 
results. Noble [ 57 ] in the same manner as Mach, was able 
to test subjects' perception of the vertical for specific 
force vectors making an anqle of from 6 to 40 degrees with 
respect to the vertical. Noble concluded; "The empirical 
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function corroborates Mach's original hypothesis; vis, , that 
Ss adjust in accordance with the resultant of centrifugal 
and gravitational forces." will find that this hypothesis 

is not generally true except for type M stimuli, 

Clark and Graybiel [ 11 ] conducted an experiment similar 
to that of Noble’s except that the centrifuge was slowly 
accelerated up to its maximum rate of rotation while the 
subject continuously tracked his perceived vertical. If the 
small tangential acceleration and the dynamic effects of the 
slowly changing stimulus is ignored then this experiment can 
also be classified as type N. • Clark and Graybiel described 
their data as follows: "The data also show that when radial 
acceleration is increased slowly , as in this experiment , the 
subject adjusts to the resultant force accurately, with a 
slight tendency for the perception to lag behind the stimulus 
and a tendency to set the line at a position slightly greater 
than 0. Insofar as the accuracy of the settings is concerned 
the data support the results of .. .Noble . w Other category N 
experiments were conducted by Clark and Graybiel [12, 30 ] 
with similar results. Thus, for these type N experiments, the 
general conclusion can be drawn that the perceived vertical 
is aligned approximately with the specific force vector. 

Schone [ 63] conducted a thorough set of experiments using 
a centrifuge to chafrge the magnitude of the specific force 
vector, SF. His experiments differed from those described 
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above in that he did not put his subjects in an erect 
position relative to the earth vertical but varied their 
orientation systematically with respect to SF. In doing so 
he conducted some trials in each of the three categories. 

i 

Figure 7.2 summarizes Schone’s data for the case in which 
the subject's roll angle (lateral tilt) is varied in a 1 g 
and in a 2 g environment. The experiments ih 1 g are all 
category A experiments (except for the case with zero roll 
angle which is category N)‘. For the experiments in 2 g 
the z component of specific force is given by 

SF- =-2cos (<j>) cos (0) (7.4) 

z 

where 0 = roll angle of ftead with respect to SF 

and 0 = tilt angle pf the utricular macule with 
respect to the sagital plane (approximately 25-30 
degrees, see Figure 7.1) . 

Therefore the experiment is of * category E if cos (0) > .5 
( cjj ^ 60°) , category N if cos (0) « .5 (0 = 60°) and category 
A if cos(0) < .5 (0 > 60°). Referring again to Figure 7.2 
it is clear that all category A experiments exhibit the Aubert 
effect, all category N experiments exhibit essentially no 
bias, and all category E experiments exhibit the Muller effect 
Schbne repeated the tests in one g under water to test for 
the possible influence of the tactile sense on perception of 
orientation. The results of these tests are in general 
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Figure 1.2 Perceived Tilt Angle as a Function of Actual 
Tilt Angle for 1G and 2 G Specific Force 
Environments 
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agreement with the lg curve in Figure 7.2 with the only 
difference being a slightly greater Aubert effect for 0 
greater than 90°. 

Finally Schone varied the pitch orientation (a) of the 
subject relative to the local direction of SF for different 
specific force magnitudes. Subjects were instructed to 
adjust the vertical position of a luminous spot until it was 
"on the horizon" or "until he sees it at eye level." Figure 
7.3 summarizes the data from these experiments (constant 

i 

errors have been eliminated by defining the perceived pitch 
angle when the head is erect in 1 g as zero) . 

To categorize these experimental trials it is necessary 
to calculate the z component of specific force. 

SF (a, G) = -Gcos(e-a) (7.5) 

where 0 = 25-30 degrees 

a = angle subject is pitched down (Figure 7.3) 
and G = magnitude of the specific force vector (g) 

All experimentally determined data points shown in Figure 
7.3 (G^l) have associated values of G and a which satisfy 

SF (a , G) < - cos (0) (7.6) 

and therefore the experiments which gave rise to them must 
be characterized as type E. To interpret these experimental 
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Figure 7.J Perceived Pitch Angle uu a Function of True’ 
Pitch Angle and Magnitude of Specific Force 
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results it is convenient to refer to the subject’s perceived 
direction of "down/ 1 which can be derived from his indications 
of the horizon, and to compare this with the direction of the 
specific force vector. All data points confirm to the 
following rule: the true direction of SF is always between 

the perceived direction of "down 11 and the negative z axis 
(~i ) . This is equivalent to a Mtlller effect in pitch with 
respect to -i z . In other words the deviation of the per- 
ceived down from -i will always be in the same direction but 
of greater magnitude than the actual deviation of the specific 
force from -i for a type E experiment. This effect is almost' 
nonexistent for the lg curve but becomes progressively 
greater as the inequality given in 7.6 is strengthened. The 
fact that all the curves cross at a~ 6 is to be expected since 
at this point -i , SF and the perceived direction of "down" 
will all coincide. Finally the data points for 1.9g, 1.6g, 

1.3g and lg can be extrapolated $t each value of 0 to yield 
a prediction for the perception of down at .5g and Og. In 
these cases 

SF (a,. G) > - cos (0) (7.7) 

and thus the corresponding experiments would be of type A. 

This extrapolated data conforms to a rule just opposite to 
that previously stated. The rule governing this data would be: 
the perceived direction of "down" is always between ^i_ z and SF. 



This is equivalent to an Auberfc effect in pitch with 
respect to -i . Data from Cohen [ 17] for which the head 
was erect in l.Og, i„25g, 1.5g # and l-75g is also shown 
in Figure 7.5, Cohen's data supports the progressive in- 
crease in the error in perceived pitch orientation as the 
magnitude of SF increases but does not show as large an il- 
lusion of pitch as does the data, from Scheme* 

For the experimental evidence ^reyiewed to this point we 
can infer the following correlation between stimulus type 
and the resulting perception of - the Vertical: 

Category A -* underestimation of tilt angle with 

respect to -i (Aubert phenomenon) 

2 

Category N ■* perceived vertical coincides with SF ^ 

Category E -> overestimation of tilt angle with 

respect to -i (Mlilier-like phenomenon) 

amo dj 

Data from Miller and Graybiol [ 52 ] along- with the 
experimental data described above will be compared in 
Section 7.3 to the predictions of a model (based on the above 
correlation) developed in Section 7.2. 
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7.2. Model Based on Altered Saccular Information 

The simplest mathematical model for the steady state 
perception of orientation with respect to a specific force 
field based upon information from the otoliths would be 
a model consisting of three orthogonal accelerometers whose 
outputs correctly. reflected the components of the true specif- 
ic force vector. Using this model, the predicted perception 
of the direction "down" would always coincide with the 
direction of the steady state specific force vector. The 
experimental results presented in the previous section 
demonstrate that, while such a model would adequately pre- 
dict responses in lg for cases in which the head is tilted 
less than 20 degrees from the erect position, this simple 
model is inadequate for stimuli not meeting these criteria . 

Before developing a model which accounts for perceptions 
of "down" which deviate from the steady state direction. of 
S F it is important to illustrate how the responses of afferent 
nerves associated with different morphological polarizations 
can be combined simply to give an estimate of SF. The 
morphological polarization maps, for the utriculus and 
sacculus of the squirrel monkey were illustrated in Figure 2.6 
Fernandez, Goldberg, and Abend [ 24 ] demonstrated that the 
response of a given otolith neuron was a linear function of 
the component of SF in only one direction. The direction of 
sensitivity is designated, the functional polarization vector F 
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Fernandez et aL found that thp distribution of the vectors 
F for the neurons they studied agreed closely with the 
morphological maps shown in Figure 2.6. If the functional 
polarization vector of the i— ceil were given by 


li - + Y i% + z O* 


(7.9 


where 


+ Y. 


+ Z. 


and , iy and i^ are defined as in section 7.1, then the 


steady state response of this cell to a specific force 

i 

stimulus SF would be given by 


FR. = S.(SF-F-) + SFR, (7.10) 

1 X A A 

th 

where FR i is the afferent firing rate of the i — • cell 
(ips) 

S. is the sensitivity of the x— cell (ips/g) 
and Sfr. ig the tonic firing rate of the i — cell (ips) . 

The firing rates from a group of such cells (i®l, N) can be 
combined to yield a signal which is directly related to the 
component of specific force along each of the three axes 
L x l iy' an ^ i 2 • For example, an estimate of the component 
of specific force along i would be given by 


SF 

z 


^{(FR.-SFR^F, i 2 )/S.} 


N 

l 

i®l 


(F, 




2 


(7. 11) 
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A simple model for the perception of "down" can then be 
formulated as follows: 


DOWN 


IWSFyiy+SF^ 
(SF^+SFy+SF^) 1/2 


(7.12) 


where DOWN is a. unit vector in the perceived direction 
of "down" 

A 

SF 2 is calculated from 7.11 
and SF X and SF^ are calculated from 7.11 with "z" 
replaced by *'x" and "y" respectively. 

Since this simple model implies that 


e[ DOm “ sf/|sf| ] » 0.0 


(7.13) 


where L[ • ] denotes the expectation operator, 
none of the deviations noted in section 7.1 will be pre- 
dicted. The question of how to modify 7.12 so that the 
various phenomenon associated with the perception of tilt can 
be predicted can now be answered. 

To motivate the model which will follow it should be 
noted that the classification pf stimuli developed in 
section 7.1 does not depend on the components of SF along the 
i x or iy axes. In other words, a change in either SF x or SF^ 
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will not alter the category of the stimuli and consequently 
will not change the type of response (Auberfc, Muller or no 
deviation) predicted* It thus seems reasonable to assume that 
the parameter which controls the deviations which we endeavor 
to predict is SF and not SF or SF * Using this assumption 

u A y 

A A 

we conclude that SF and SF ' should be calculated in 

x y 

accordance with equation 7.11 apd they should then be used in 

equation 7.12 without alteration. SF , on the other hand, must 

z 

be altered to predict the illusory perceptions found experimen- 
tally. While it is relatively easy to motivate the altera- 

A 

tions of SF qualitatively, the agreement of the model's 
6 

predictions with the experimental data will serve as the 
major justification for the specific alterations proposed in 
the mode h 

Figure 7.4 illustrates stimuli from the three stimulus 

categories 0 Since for category ^ stimuli w© want the model to 

predict that the direction of SF will be correctly perceived, 

A 

we should use the value for SF calculated- from 7.11* If 

z 

A 

the value substitued for SF in equation 7.1’2 is designated 

Z 

SF Z then our first rule is that 

1* If SF = -cos (9) (category N) (7.14) 

Z 

/*W A 

then set SF z = SF z 

For a category A stimulus we should predict a perceived 

direction of "down" between SF and -i . This will be 

— — z 

accomplished if SF is decreased. Therefore 

z 




( 7 o 15) 


2= If SF > -cos(9) (category A) 

Then SF < SF 
z z 

For a category E stimulus we must predict a perceived 
direction of "down *' which deviates from bv more than 

/V 

SF. This will be achieved if SF„ is increased (made less 

z 

negative) <, Therefore 

3 0 If SF r < -cosO) (category E) 

2 V ’ 

Then set SF > SF 

Z 2 


Figure 7.5 presents the exact alteration of SF^ which 
is necessary to fit the experimental data. Figure 7.6 shows 
the resultant model for perceived orientation. The only 
difference between this model and the simple model is the 
nonlinear alteration of SF . Xn this model the accelerometers 
are presumed to have unity gain, Since most of the otolith 
neurons with functional polarization vectors which lie close 
to the xy plane are utricular in origin and since most of 
those which have significant components along i originate 

' mj Z 

from the sacculus the accelerometers in Figure 7.6 along 
— x an< ^ —y are Isheled "utricle" and the accelerometer with 
sensitivity along i is labeled "saccule." The coordinate 
transformation is inserted to take cognizance of the fact that 
our perception of "down" is usually related to the head’s 
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principle axes and not to i , i , i . The coordinate trans- 
formation and normalization is given by: 


DOWN xhd 


cos ( 0 ) 

0 

- sin(0) 


SF x 

DOWN yhd 


0 

1 

0 


A' 

SF y 

down 2HD 


sin (0) 

0 

cos (9 ) 



. * 




. 


. 


(7.17) 


2 — 2 ' 

SF z + SFy + SF z 


1/2 


where i v , m = i (head erect pitch axis? pitch forward 
~ YHD “ y positive) 

izw ~ yaw axis head erect {yaw left positive) 

ixHD * iyHD x izHD ( head erect roll axis; roll 
* ~ right positive) 

and Q -tilt angle of the average utricular plane with 
respect to sagital plane (25° will be used in 
section 7.3). 



i 



7. 3o Model Predictions in Various 
Environments D _ ~ 


rce 



If the four conditions discussed in section 7.1 are pre- 
sumed to be met then a given experimental trial can be fully 
specified by the specific force vector? SF, SF will be re- 
ferred to the standard head fixed coordinate system (see 
Figure 7 .7) 0 The orientation ©f §F can moat easily be pic- 
tured if the angle between ths median plane and a vertical 
plane containing SF is given (fl-gp) and the angle- between 
SF and -i^is given C^ SF ) ° Once the orients tioro of §F is 
determined the only parameter which must be known is the 
magnitude of SF which will be designated by 

The problem of how to interpret the predicted perception 
of “’down" is somewhat more difficult. The implicit assump- 
tion made by most experimenters has b eon that the perceived 
direction of "down" will always be in a plan© which includes 
both i and SF, In fact from the evidence given in the last 
two sections it seems more likely that th© vector DOW? lies 
in a plane which includes i and SF. Typically*, in a lateral 
tilt experiment (e.g.? 90° tilt right in Ig) a subject 
would be asked to adjust the orientation of a luminous line 
in the frontal plane until it appeared vertical but no test 
would be given to see whether in this orientation he might 
not feel that he is oriented slightly face down or face up. 
Experimental data from Nelson [ 56 ] demonstrated that when a 
subject attempts to position himeelf to a 9.0® roll right 
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Figure 7.7 Orientation of SF With Respect to 
Head Axes 



Figure 7.8 Orientation of DOWiSi with Respect to 
head Axes 



orientation underwater he actually ends up positioning himself 
slightly ( a S°) face up» Sine© this deviation was £ound to 
bo statistically significant w© can assume that in an experi- 
ment in which the subject is just rolled 90° he must sense 
that his orientation is slightly face down* This deviation is 
consistent with the model shown in Figure 7,6, The problem 
which naturally arises is that if DOWN is specified in the 
same way as Sj? (by giving values for and jj ) then it 
cannot be readily compared to the experimental data* Figure 
7„8 shows how this comparison is simplified. For an experi- 
ment in which the subject is asked to indicate his perception 

of the vertical by adjusting the orientation of a line in his 

frontal plane,, the data is compared to the angle (see 
Figure 7 0 8 ) which is the angle h© must roll in order to put his 
perception of "down” into his median plan®* For an experiment 
in which the subject is asked to indicat® his perceived pitch 
angle , the data is compared, to the angle which is the 
angle he must pitch backward to {a ring his perception of "down" 
into the iyhD—ZHD P-^- ane ° In any experiment in which SF 
is in the median plane this distinction is unnecessary since 
R d will equal zero. 

Figure 7.9 summarizes the model * s prediction for the 
perception of lateral tilt (hg F = 90°) in a Ig and a 2g 

environment. The value of is plotted as a function of 
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Figure 7 . ‘J Model .Predictions for Perceived Tilt AngU 
as it Function of Actual Tilt Angie in 1C 
and 2G Specific Force Environments 
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Vi for G„ „ * 1 and G_„ ® 2. The data from Schone [ 63 3 
£> r i* i? r 

which was given in Figure 7.2 is repeated along with data 
from Miller and Graybiel I 52 3 • The small arrows indicate 

. > /s 

how these predictions would change if the alteration of SF^ 
shown in Figure 7.5 for type E stimuli is increased. 

Figure 7.10 summarizes the model's prediction 
for the perception of pitch (^gp =0°) in Ig and 

2g environments. The value of P Q is plotted as a function of 
]i„„ for G__ = 1 and 2. Most of the experimental data avail- 

oJ 1 bf 

able for the perception of pitch orientation was gathered 
for relatively small pitch angles (P gF < 40°). Figure 7.11 
shows the model's predictions of as a function of Vi gF 
in the region from 0 to 40 degrees for G gF = 0g> . 5g, 
l.Og, 1.3g f 1 . 6g and 1.9g. These curves should be compared 
to the data for Schone and Cohen (Figure 7.3). The predicted 
deviations of the perception from the trues pitch angle varies 
as a function of l J gF and G gF in exactly the same way as does 
Schone’ s data except that the predicted deviations are only 


about 60S as great for the cases in which G gF > lg. On the 
other hand, the model accounts for 90% of the pitch percep- 


tion exhibited by Cohen’s data. These discrepancies can be 

remedied if the alteration of SF is increased for SF < -cos ( u ) 

z z 

(category E stimulus) as shown ir> Figure 7.12. A slope of -0.4 


is needed to fit Schone' s data while a slope of -0.1 is 


sufficient to fit Cohen’s data. The small arrows perpendicular 
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Figure? 7.10 Model Predictions for Perceived Pitch Anglo 
as a Function of Vruc Pitch Angle in 1G and 
2C, Specific, Force Environment si 


i 


PERCEIVED PITCH ANGLE (DEGREES) 
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to the curves in Figures 7.10 and 7.11 Illustrate the direc- 

A 

tion in which the predictions move as this alteration of SF z 

is increased. Finally , the model predictions were generated 

using a value of 25 degrees for 9 (the angle between the 

average plane of the utricles and the i axis ) while the 

XnD 

data from Schbne would be best fit if an angle of approximately 
29 degrees were used. Since 29 degrees is within the 25-30 
degree range normally associated with the orientation of the 
utricles this does not constitute a discrepancy . 

All category IS! stimuli can be su aimed up by stating that 
the model correctly predicts that the perception of the verti- 
cal will be aligned with the specific force vector with essen- 
tially no error. The experimental justification for this was 
discussed at length in Section 7.2 [11,12/30,45,46,47,57]. 

while it might bo noted that all model predictions and 
all experimental results which have been given are for cases 
in which the specific force vector was in the median plane 
(pure pitch) or in the frontal plane (pure lateral tilt) the 
model is capable of generating predictions for an arbitrary 
constant specific force vector. 
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7.4. Summary 

The illusions associated with the perception of static 
tilt in various specific force environments have been re- 
viewed and then classified in such a way that a simple per- 
ceptual model could be developed to account for the 
available experimental data. The fundamental conclusion to 
be drawn from this model is that these illusions can be 
accounted for by a simple nonlinear transformation of the 
information (mainly from the saccule) related to the component 
of specific force perpendicular to the average plane of the 
utricles. • 



CHAPTER VXII 


INTEGRATION OF SEMICIRCULAR CANAL MID OTOLITH INFORMATION FOR 

MULTISENS QRY STIMULI 


In this chapter we consider the class of stimuli speci- 
fically excluded in. Chapter Three? namely those stimuli which 
simultaneously excite the semicircular canals and the otoliths 
in generalizing the stimulus class to include any combination 
of rotational acceleration and translational acceleration in 
three axes, a number of significant new problems arise*, These 
problems and the philosophical approach taken to deal with 
them are discussed in section 8.1. A mathematical model of 
the perception of dynamic orientation which results from the 
combined effect of arbitrary angular and translational accel- 
erations is developed in section 8.2. The model * s qualitative 
predictions for several stimuli are also discussed. Quanti- 
tative predictions from a computer simulation of the model 
are presented for several stimuli in section ‘8*3 along with 
the model 5 s frequency response for small pitch and roll angle 
oscillations . 
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8 . 1 Discussion of Modelling Problems and Philosophy 

Before discussing any of the problems associated with 
the integration of sensory information from the semicircular 
canals and otoliths it is important to clarify what the im- 
portant perceptual outputs of the model should be. Certainly 
we would be interested in estimating the following quantities: 

1. Orientation of the head with respect to the 
gravitational vertical 

2. Rate of rotation of the head about its three , 

principle axes 

3. The translational acceleration of the head with 
respect to its three principle axes 

and 4. Additional quantities which are derived from 
the preceeding (e.g. azimuth/ translational 
velocities and translational positions). 

The most important of these is the determination of 
orientation with respect to the vertical. Strictly speaking, 
there is no way of using information which is derived only 
from the otoliths to determine , the direction of the gravita- 
tional vertical if there is no a priori information regarding 
the expected variations in orientation or translational 
acceleration. The principle of equivalence in general rela- 
tivity precludes such' a separation based purely on measurements 
taken from linear accelerometers.' In the language of modern 
control systems theory any system comprised only of linear 
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accelerometers is unobservable • How then are we capable of 
distinguishing a change in orientation with respect to the 
gravitational vertical from a change in acceleration? The 
answer to this question has two parts ► Firsts we are not 
restricted to the use of linear accelerometers (otoliths) 
since we also have angular velocity transducers (the semi- 
circular canals) which indicate with reasonable accuracy 
the rate of change of the head 5 s orientation for rotational 
rates in the frequency range from 0.1 rad/ sec to 10 rad/sec. 
Roughly speaking , for changes in the direction of specific 
force which occur in this frequency range (as determined 
from otolith information) the distinction between a change 
in oreintation with respect to the gravitational vertical 
and a transla tional acceleration (or some combination of the 
two) can be made by noting the output of the semicircular 
canals* As the frequency of the variations in the direction 
of the specific force vector decreases below <, 1 rad/sec , 
information from the canals becomes less and less useful. 

In fact as the frequency of these variations approaches zero 
the system is incapable of determining the gravitational 
vertical. The second part of our answer therefore is that 
for lower frequency variations the system cannot concern 
itself with the true gravitational vertical but must be 
content to estimate an "effective graviational vertical" which 
can serve as the practical reference for man's normal 
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activities. The phenomenon of associating the gravitational 
vertical with the perceived direction of specific force for 
very low frequency (essentially static) stimuli was discussed 
in depth in Chapter seven. 

Once the direction of the gravitational vertical is 
estimated, the other perceptual quantities can be derived. 

The sensation of rotation about an axis parallel to the per- 
ceived gravitational vertical : ) will reflect exclusively 
the properties of the semicircular canals described in Chapter 
Three . 

The perception of rotation about an axis perpendicular 
to the perceived vertical (w^) should reflect the information 
available from the canals and th§ otoliths. Since the otoliths 
are capable of sensing a constant change in orientation with 
respect to the gravitational vertical. the perception of con- 
stant rotation about a horizontal axis • should persist indefin- 
itely. Henson and Bodin (Ref. 4 ) and Guedry (Ref. 34 ) confirm 
that the perception of rotation does indeed persist for prolonged 
rotations about a horizontal cephalocaudal axis. 

The estimate of translational acceleration is essentially 
determined once the direction of gravity is estimated since 

A ~ G - SF (8.1) 

where A = translational acceleration (g’s) 

G = gravitational vector (normally 1 g) 

SF = net specific force vector (g's) 


and 
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The only change needed in equation 8,1 is the replacement of 
each term by its estimate (e»g 0 A by A p etc 0 )„ To maintain 
the notation used in chapter Seven the estimate of <3 will be 

/N 

denoted by DOWN since this is more descriptive of its per- 
ceptual meaning.- The remaining perceptual quanitities (azi- 
muth, translational velocity and translational position) are 
obtained by integration as follows? 


?(t) 


*<t 0 ) + 



( 8 . 2 ) 


v(t) = v(t Q ) + t / A(n)dn 


x(t) = x(t Q ) + t / v(n)dn 


( 8 . 3 ) 


(8,4) 


where V is the angle between the • projection of the head’s roll 
axis (i. XHD ) in the earth 8 s horizontal plane and some 
fixed direction in that plane (e.g. a vector pointing 
toward true north) . 

A 

is the perception of rotation about an axis parallel 
/\ 

to DOWN 


A 

V is the perception of linear velocity 

A 

and x is the perception of spatial position. 

In ally the model should be capable of predicting 15 

• , A A 

quantities (3 associated with w, 3 associated with A# 3 asso- 

. . a A 

ciatcd with V, 3 associated with X,. 2 associated with the 

A A 

direction of DOWN , and 1 associated with ¥) . Of these 15, 

A 

the 2 associated with the direction of DOWN are by far the 
most difficult to model quantitatively and for this reason 
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Ui<- inorlo.1 developed in this chapter is called the "down" 
estimator. Equations 8.1, 8.2, ,8.3 and 8.4 determine the 

A * A A- ^ 

quantities A, 'F , V, and X as a function of DOWN and w )( 

and therefore the only estimators which are left to be 

^ A ^ 

modelled are those which generate DOWN and to (since cu ( ^ is 

A ^ 

the component of to parallel to DOWN ) . 

Before considering these estimators (for DOWN and w) in 
detail, several problems require consideration. The first of 
these is the problem of reconciling what may seem to be con- 
tradictory information from the canals and otoliths. Three 
examples can be cited for which there exists corresponding 
data.. The first of these involves an abrupt change in the dir- 
ection of the specific force vector relative to the head ("rot- 
ation information" from the otoliths) without any correspond- 
ing indication of rotation from the canals (e.g. aircraft 
catapult launch Ref. 18 , or a change in the direction of 
specific force due to rotation on a centrifuge Ref. 28 ) . 

A second example of such a conflict would arise in the case 
of a constant rotation about a horizontal axis which would 
lead to a continuously rotating specific force vector but 
a zero steady state output from the canals (e.g. a barbecue- 
spit experiment Ref. 4 , 5 , 34 ) . Finally, situations may 

arise in which the canals indicate an abrupt rotation about 
a horizontal axis but the otoliths indicate . no change in the 
direction of specific force (e.g. a coordinated aircraft turn 
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or the abrupt cessation of rotation in a barbecue-spit experi- 
ment Ref 4 B 5 * 34 } ° Since one or more examples from 
each of these categories will be treated in detail in the 
remaining sections of this chapter it is unnecessary at this 
point to give a full accounting of the perceptual responses 
except to say that the perception of the vertical for these 
stimuli is most strongly associated withs 

1. The low frequency portion of the ‘'rotation infor- 
mation from the otoliths 

plus 2. that part of the canal information which is 

consistent with the high frequency portion of the 
"rotation information' 1 from the otoliths. 

Since the rate of movement of the perceived vertical may 
not be consistent with the estimate of rotation based only 
upon canal, information the question arises whether the per- 


/V 


caption of rotation reflects the movement of DOWN or canal 

informal, i on or a combination of the two. If the time histories 

A a -a 

of DOWN ami d'j (the component of fu perpendicular to DOWN ) were 

to be consistent then in the situation in which the direction 

^ ✓A 

of DOWN is constant it must follow that = 0_. The experi- 
mental evidence' (Ref. 4 , 5 , 34 ) does not consistently 

support this conclusion and thus DOWN and may not be in 
agreement. Although such a contradictory sensation (of rot- 
ating but not changing one's position) seems difficult to 
imagine, it :i.s also found in cases in which otolithic and 
visual information conflict (Ref. 21 ) and during caloric 
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testing. The fact that these sensations are contradictory 
also complicates interpretation of some of the experimental 
data. For example, if an experimenter asked a subject if he 
felt himself rotating the subject could answer either "yes” 
or "no" (in fact an answer of yes and no would be more 
appropriate ! ] . 

A second problem arises in the case of stimuli which are 
predictable, usually because the subject is thoroughly familiar 
with the stimulus from past experience and is able to recog- 
nize the underlying stimulus pattern. The phenomena associ- 
ated with such a situation are significantly different then 
those which we are attempting to model here since they in- 
volve the complex problems of pattern recognition. Furthermore, 
it is very likely that the processes involved in recognition 
are strongly dependent on the simplicity of the stimulus, the 
subject's past exposure and many other factors which would 
make an accurate prediction of the perception of dynamic orien- 
tation extremely difficult. For these reasons the stimulus 
class for which we are at tempting to model the perceptual 
responses will be assumed to be unpredictable. 

Finally, the information upon which the "down” estimator 
bases its estimate must be considered. Although the informa- 
tion from two sets of semicircular canals and otoliths is 
available to the brain it is unnecessary to waste computation 
time performing a dual set of sensory simulations.. For this 
reason, the model simulates cyclopian sensors located near the 
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center of the skull. One case in which this simplification 
could make a difference is one in which, a very rapid rotation 
of the head is made about either the yaw or roll axis* in 
which case the otoliths on either side of the head would nor- 
mally be exposed to slightly different specific force vectors. 
Since no perceptual effect has been ascribed to this difference 
nothing is lost by eliminating it* Figure 8*1 illustrates 
the orientation of the three canals which are aligned with and 
sensitive to rotations about the three axes i * i vn and i , 

These axes are equivalent to the axes i * t and i„ illustrated 

—x — y ■— z 

in Figure 7,1 except they are rotated 4-45° about i = i 7pB The 
canal labeled LH has an afferent response equivalent to that 
of the left horizontal canal and opposite to that of the right 
horizontal canal. The canal labeled RS has a response equi- 
valent to that of the right superior canal and opposite to 
that of the left posterior canal. Finally the canal labeled 
1.53 has a response equivalent to that of the left superior canal 
and opposite to that of the right posterior canal. Each canal 
and the optimal, estimator for the rotational rate about its 
associated axis is modelled by the dynamic model and Kalman 
filter developed in section 3.1. The otolith models for the 
utricle and saccule along with their associated processors 
which were developed in section 3.2 are used to generate the 
estimate of specific force. One utricular sensor is aligned 
with its major axis of sensitivity along i v = i and one along 
ivo = — Y ^ seo Figure 7.1). The saccular sensor is aligned with 
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Figure 8.1 Orientation and Sensitive Axes of 

Cyclopian Semicircular Canal System 
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i - i„ = i„ n (see Figures 7*1 and 8*1) * 

Figure 8o2 summarises the information available to the 
"down" estimator* The alteration of saccular information shown 
in Figure 8*2 is identical to that developed in Chapter Seven 
and its presence insures that the steady state performance of 
the "down" estimator for static tilts will be consistent with 
the predictions made in that chapter* The estimates of rota- 
tion based upon canal information and specific force based 
upon otolith information are transformed from sensor to head 
coordinates before being used by the "down 08 estimator since 
the principle head axes are the most natural coordinates to 
which to refer our concious perceptions of dynamic orientation* 
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8.2 DOWN Estimator 


After numerous algorithms were developed in an attempt 
to successfully produce an estimator with the desired quail- 
fications, one was found which fulfills all of the require- 
ments discussed in section 8 . 1 and which yields very reason- 
able quantitative results (see section 8o3). The discrete 
"down" estimator is illustrated in figures 8.3 and 8 0 4 „ The 
information available to the “down* 1 estimator at the begin- 

A 

ning of each update is the old estimate of down, DOWN ( t-At) 

A 

and the new estimates based upon canal information (^^(t)) 

■ A 

and upon otolith information (SF^Ct)). Figure 8-. 3 illus- 
trates the calculation of the updated perception of down, 

A 

DOWN (t) and Figure 8.4 illustrates the updated perception of 

A 

rotation, co (t) « Each element of the model is labeled with 
a letter from A to L for easy reference and will be discus- 
sed in alphabetical order.. 

The first element, labeled A and marked with an X, 
represents the following computational procedure? Produce 

A A 

a vector which is in the direction SF tTy> (t-e)X SF U ~ (t+e) 

and which has magnitude equal to the angular rate of change 
of the direction of SF^ p at time t. In the computer simu- 

A 

lations carried out in section 8.3, w ur . was calculated each 

1 — fiD 


second (t = 1,2,3, 


•) and SF^p was calculated on the 


1 c Am . 

half second (t = j j j . . . ) so e = 2 ~ = *5 seconds « 
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Figure 8,3 DOhH Estimator 
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Figure 8.4 
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represents the information available from the otoliths con- 
cerning the rate of rotation of the head i£_ it were assumed 
that SF was fixed in space. 

The low pass filter, labeled B, performs the function 
of separating out the low frequency component of which 

is assumed to arise from the change in the body’s orientation 
with respect to the gravitational vertical. The output sig- 
nal ajk is intended to fill in the low frequency information 
—Sr 

A 

missing from the canal signal uj_ hd for rotations about a 

horizontal axis, is the high frequency component of 

—or 

ui^and typically arises from both transient linear acceler- 

— or 

ations and abrupt changes in the head ' s orientation with 
respect to the gravitational vertical. The best time con- 
stant for the low pass filter was found to be approximately 
35 seconds. 

The transformation labeled C produces a rotation vector 
from Wgp as follows: 


^OTO 


Component of which is perpendicular to 

the plane of SF and DOWN ( t-At) 


It may seem odd at first that this transformation allows 

A A 

rotations which would by themselves move DOWN away from SF . 
The reason for this is that such rotations are necessary to 
cancel the canal signals which arise when prolonged rotations 
are suddenly stopped. It is this mechanism which helps to 
predict the stabilization of the perception of orientation 
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when prolonged rotations about a horizontal axis are abruptly 
terminated as was found experimentally by Benson and Bodin 
(Ref o 4,5) and Guedry (Ref„ 34 ) . In all the simulations 
carried out no case has been encountered in which Rg CC (which 
will be discussed next) did not cancel completely any R QT0 

/N A 

which would move DOWN away from SF. If such a case occurred 
it would appear reasonable to decrease the magnitude of Rqijiq 
until the net effect of = 5 qT0 ~ —SCC wou ^ d ke to 


mize the misalignment of DOWN and SF after rotation . The 
combined effect of elements A, B and C in figure 8.3 is to 
produce a* rotation vector Rq T 0 from the current and past 
estimates of? SF which represents the low frequency rotational 
rate information due to the otoliths » 

We now turn our attention to the information available 
from the canals. The rotational information from the semi- 
circular canals must be consistent with the high frequency 
sensations arising from the otoliths (represented by Wg p ) 
if it is to be used to update the sensation of orientation 

A 

with respect to the vertical. The portion of co^ (t) which is 
consistent with w sp is denoted by and is calculated by 
the following procedure: 


1 ) 

2 ) 

3 ) 


Calculate the component of which is paral- 


H 


lei to <*>__. Call this component C. 

— of — 


H 


If C is in a direction opposite to w_g p then set 

£ 

H 

if C is in the same direction as 


then set 
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C 


~ W SF if l£l > 


^SF 




= ~C 


if |C 


i4f 


( 8 . 6 ) 


A H 

The portion of to 4ir ^(t) which is inconsistent with is 

— rlL’ -~br 


denoted by and is given by 


Sc = *VD lt) - ^ 


(8.7) 


The reason that -w HD (t) (instead of o^t)) was compared 
H 

with co gF is that for a positive perception of rotation the 
corresponding rotation of the g vector would be negative. 

While experimental evidence clearly indicates that the 
effect of on the perception of orientation is minimal it 
is not clear that it has no effect in the very short term 
(<1 sec) , For this reason is passed through a high pass 
filter (E) of the form t s/(;s+ 1) where t < 1 sec. For the 
catapult launch simulation described in section 8.3.3 t 
could be no higher than .25 seconds to retain reasonable 
results. A value of 1=0 would not be inconsistent with 
any available experimental evidence. 

The rotation vector due to canal information is denoted 

. C 

by R . R ,,, is computed by taking the sum of and the 

result of filtering and then eliminating the component 

which is parallel to the last estimate of down (since this 

A 

component is ineffective in changing the direction of DOWN 
relative to the head) . 
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is then computed by subtracting Rg CC from . 

represents the estimate of the rotational rate of the 
outside world around an axis perpendicular to the last esti- 
mated direction of down for the purposes of updating that 

A 

estimate*. The transformation labeled G updates DOWN ( t- At ) 
using . The output of G is denoted by D ' (t) and 


satisfies 


1) DOWN (t-At) x D ' it) is in the same direction as ^tq T 

and ^ 

2) the angle between DOWN (t~At) and D° (t) is given by 

I Spot I At 

( 8 , 8 ) 

Therefore, if ISjq^I = 30 degrees/sec. and At = 0.5 seconds 
DOWN will be rotated about 5p QT by 15 degrees. 

D # (t) would normally be considered the new estimate of 
"down'* except that because it is generated through the in- 
tegration of rate information it is bound to accumulate errors 
which must be eliminated if permanent discrepancies are to 
be avoided. This is accomplished through a slow reduction 

A 

of any discrepancy in direction between D ! and SF^ (elements 
H and I). The time constant, i p , is quite large but was 
found to be a weak function of the magnitude of the specific 

A 

force vector (as | SF | increases, Tp decreases and D‘ moves 
toward SF more rapidly) . Tp is given by 
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60 

ITT 17 * 


(8.9) 


The net effect of H and I is that in the steady state 
the subject adopts the estimated specific force vector, 
based on otolith information, as the correct direction for 

A 

DOWN. This insures that the steady state response of the 
model will exhibit the perceptual errors modelled in chapter 
Seven . 


The resulting estimate of DOWN (t) represents the model's 
prediction for the subject's perception of the direction of 
the gravitational force vector with respect to his head. 

This estimate is then used at time t + At to generate a new 
estimate . 

The model for predicting the perceived rate of bodily 


rotation is shown in Figure 8.4. to^ <t) is found simply by 

A A 

taking the component of <o HD (t) parallel to DOWN (t) . wg is 
defined to be the bodily rate of rotation which would be 

A 

consistent with the rate of change of the direction of DOWN . 
The transformation K is similar to that at A in figure 8.3 
except for a minus sign. <o^(t) is formed by: 

1) calculating the difference between the component 

A A 

of (t) which is perpendicular to DOWN and tog 

2) passing this difference through a high pass filter 


(L) 

and tnen 

3) adding the resulting output to top- . 



This arrangement accepts the relatively high frequency changes 
in rotation rate indicated by the semicircular canal system 

A 

while deferring to the rate of rotation consistent with DOWN 
for lower frequency changes. Data from Benson and Bodin 
(Ref. 4,5) indicates that a filter of the form ^/(tS-J-I) 
with t - 0-5 seconds should be sufficient (if t « 0 then 

A ^ A 

tu (t) = wg and Would be consistent with DOWN ) 0 The total 

A A 

sense of rotation, o(t), is given by the sum of oj . (t) and 

— 

(t) . 

This completes the component by component review of the 
model. Before describing the quantitative results which were 
produced by .computer simulation* several examples of qualita- 
tive predictions will be giv®n<, 

First consider a standard rate aircraft turn which is 
abruptly stopped by rolling out of the turn rapidly into level 
flight. Just before the rollput the subject will perceive 
himself to have zero roil angle with respect to the earth 
vertical and a slightly pitched back orientation due to the 
slightly increased g force in t^he turn (elevator illusion) . 

In addition he will have no sense of rotation since the canal 
response to the rotation of the aircraft has long since 
decayed to zero and tug = 0 . During the roll out the specific 
force vector will remain aligned with the yaw axis of the 

A 

body and diminish in intensity to lg. SF will therefore 
slightly diminish in intensity and will pitch about 1 or 2 



213 


degrees (to eliminate the slight pitched back, sensation) . 

/■N 

Since the direction of SF remains practically constant the 
otolith pathway to Rq T 0 can be considered inactive. Since 


H 


a> gF will also equal zero all of will be considered 


-HD 


inconsistent (both that part of a> generated by the rolling 

—HD 

out rotation and that due to the after sensation of stopping 

A 

the aircraft’s turn rate). Consequently all of is 

passed through the high pass filter and is quickly reduced 
to zero. Therefore, for rough calculations R = 0 and 

wLL 

A 

except for the elimination of the elevator illusion DOWN 
will remain essentially unchanged and the subject should 
sense that he is erect. 

A 

Since DOWN is essentially unchanged ^ in figure 8.4 

A 

is approximately zero. The component of which arose 

from the roll out motion of the aircraft is 1 to DOWN and 
will therefore be assigned to Since « 0 _, is set 

equal to wj and is high pass filtered with a time constant 

A 

less than 5 seconds, oj^(t) equals the output of this filter > 
(since w D = 0) which merely implies that the rolling sensa- 
tion is shorter lived (due to the high pass filter) than it 
would have been if the otolith information had not contra- 
dieted it. The component of which arose from the air- 

craft stopping its rate of turn will be in the opposite 
direction to the original turn and will be essentially paral- 

f A 

lei to the direction of DOWN . Therefore this component of 



o> will become w ^ (t) and will not be diminished., The 
sensations described above arts consistent with the illusions 
known to be associated with aircraft flight, Circumstances 
which could interfere with these illusions are the following 

1) A passenger with extensive flying experience who 
expected the turn or roll out might be capable of 
interpreting the sensations correctly . 

2) The pilot who initiated the roll out would certainly 
have little inclination towards illusions, 

or 3) Any visual Information would affect the predicted 
perception since the model presumes that there are 
no visual cues. 

A second example is that of a step in lateral acceler- 
ation of lg. Initially the subject correctly perceives him- 
self to be in an erect position in lg. Since the subject is 
never rotated during the experiment the canals are not stimu 

A 

lated and ~ 0. Referring to figure 8*4 we can conclude 

that ; 

A A -rC 

u(t) = OK (t) « U - — (t) (8.10) 

1 tS+1 

The only active pathway in figure 8.3 is that for the 

A 

information from the otoliths. S? will move very rapidly 
toward SF and then stop which will induce a rapid rise in 
u !cjp followed quickly by a rapid decay to zero, w^p will 
rise quickly during the period in which is large and 

will then slowly decay to zero. Since is perpendicular 
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to both DOWN (t- At) and SF it will pass through C and ~ 

T A 

oi . Finally DOWN will move toward SF at a rate proportional 
— - ■ 

to the magnitude of (actually a little faster since the 

lower pathway in figure 8.3 will help somewhat in moving 

A A. 

DOWN toward SF) . Figure 8.5 shows a rough sketch of the 
approximate time course of these signals. Note that since 
this stimulus is of type N (see chapter seven) the steady 
state perception of down should align itself with the true 
specific force with no error. 

The last case to be considered before presenting quanti- 
tative results is the phenomenon associated with the experi- 
ments of Benson and Bodin (Ref. 4,5) and Guedry (Ref. 34 ) . 
For a steady state rotation of o about a horizontal axis: 


and 


A 


—HD 

-*■ 

0 

^SCC 

-+■ 

0 

L 

-SF 

-+ 

-SF " 

5oto 

/\ 

ss 

5.SL 

DOWN 


SF = 


( 8 . 11 ) 


Each of these can easily be understood by reference to 
Figure 8.3 except possibly the last relation. It is clear 

that DOWN will approach SF if it is understood that the rate 

/\ 

of rotation of DOWN (a) p ) will eventually match that of SF 
since approaches the true rotation rate and any constant 

discrepancies (phase lags) will be eliminated by the lower 



Degrees, Degrees/Second 
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Figure 6.5 Approximate Time Course of Model Parameters 

and Response to 1G Step in Lateral Acceleration 
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pathway. Consequently the subject’s steady state sensation 
of rotation during the period of rotation should correctly 

A 

reflect the true rate of rotation (w » w ) . 

Immediately after the rotation stops we can predict that 
Cl 1 ' igure 8 . 3) : 


w will quickly ■+• -o>_ and then decays to zero 

(this is the typical velocity step response of 
the canals) 

0 

and then decays to zero 
-w and then decays to zero 
-w and then decays to zero 
Rq T 0 wS - 11 remain at and then cjecays to zero 

Btot = 5qto ~ —SCC ^ickly ■'* — 

and furthermore (figure 8.4): 

u)_ D will quickly ■* 0 

/X /X /V . 

w f - w, / (since aj un _l DOWN) and w „ =0 

-~X ~nL) — nO — — J| 


Wgp will quickly 

uj_ S P will quickly 
C 

will quickly 
i^ cc will quickly 


and 


tS 


^(t> = = tS+I s-j. 


( 8 . 12 ) 


(8.13) 


U) 


Therefore the model predicts- that while a subject 
should perceive that his position with respect to the verti- 
cal is not changing after the rotation ceases he may have 
(depending on the value of t chosen) a brief sensation of 
rotation opposite to the original rotation. Benson and 



Boden (Ref, 4,5) had some subjects who reported a brief 
sensation of rotation and some who didn’t. Whether this 

A 

discrepancy in reporting is due to the conflict between DOWN 
and w or due to different subjects having different values 

of t is unclear. That subjects perceive themselves to have 

A 

a constant orientation relative to the vertical ( DOWN constant) 
is not in question, Benson and Sodin report " « . . that they 
(the subjects) were quite aware that the stretcher had 
stopped and of its position relative to the gravitational 
vertical,,," Similar stimuli and reports of subjective 
responses are described in Ref, 34 . 
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8.3 Quantitative Model Predictions 


The model developed in this chapter has been computer- 
ized so that quantitative predictions can be made for arbi- 
trary stimulus combinations. The programs were written in 
Fortran XV and they include ali functions shown in figures 
8.2, 8.3 and 8.4. Although the model could be implemented 
with any update interval. At * 1 sec, was chosen as a reason- 
able compromise between computational efficiency and simula- 
tion bandwidth. One update interval takes approximately .08 
seconds of central processor time when utilizing an IBM 370- 
16 5 computer*. 

8.3.1 Dynamic Elevator Illusion 

In chapter Seven the elevator illusion was discussed 
and a model which correctly predicts its occurrence and magni- 
tude was developed. The transition from head erect in lg 
to the perception of backward tilt with the head erect in 
1.75 g was used as a test of the dynamic model developed in 
this chapter. The stimulus input to the model consists of 
a step in upward acceleration of .75g after the model was 
stabilized with head erect in lg. No rotation stimulus was 
used. Figure 8.6 shows the time course of the predicted 
pitch sensation which resulted. Superimposed on the model’s 
prediction is the data from Cohen (Ref. 17 ) in which subjects 




220 



15 30 45 60 75 

Time (seconds) 


Figure 8=6 Dynamic Elevator Illusion (1.75G) 
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were given essentially the same stimulus except that the 
acceleration was produced by a centrifuge, Cohen’s subjects 
perceived a maximum change in pitch orientation of approxi- 
mately -19°. The discrepancy in the magnitude of the steady 
state illusion is discussed in Chapter Seven. 

8.3.2 Rotation to Lateral Tilt of 5 Degrees 

Experiments were conducted in the Man Vehicle Laboratory 
by Tang (unpublished) in which subjects who were originally 
erect were tilted 5° laterally over a period of 5 seconds 
and then held at 5° for several minutes. During the entire 
period of the stimulus the subjects were instructed to adjust 
a line to their perceived vertical while being deprived of 
visual cues. The orientation of the subjects' head was con- 
trolled by using a bite board. The time course of the lateral 
tilt was recorded and used as the stimulus input for the model 
shown in Figures 8.2 and 8.3, Figure 8.7 shows the stimulus 
and the resulting prediction of the perception of lateral 
tilt. The records of subject responses were used to find the 
average peak perception of tilt and the time of its occurrence. 
In addition the average perception of tilt after one minute 
was calculated. These experimental data points are shown in 
the figure. 

8.3.3 Catapult Launch 



Tilt Angle (Degrees) 
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Figure 8.7 Perception of Lateral Tilt 
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Cohen et. al . (Ref. 18 } used a centrifuge to simulate 
the accelerations encountered in a typical aircraft catapult 
launching- The average acceleration profile used by Cohen 
is shown in Figure 8.8 along with an actual catapult launch 
acceleration profile. Figure 8.9 illustrates the manner in 
which the acceleration was generated on the centrifuge. The 
following acceleration profile was used in the simulation of 
the "down" estimator: 


A XHD * 3.8SIN( nt/3.2)g 

* o g 


t < 3.2 seconds 
t >3.2 seconds 


(8.14) 


The rotation profile used in the simulation is given by: 


0) 


2j_ 

ZHD 3 2 (1. - COS ( 2 nt/ 3.2)) deg/ sec t < 3.2 
= 0 deg/sec t > 3.2 


(8.15) 


Figure 8.10 illustrates the movement of DOWN in response to 
this stimulus. In addition to the pitch sensation for which 
Cohen et. al. tested, the model predicts a possible rolling 
sensation. If this rolling sensation is truly absent then 
the time constant in the high pass filter (element E of 
figure 8.3) should be reduced to zero. If the sensation of 
rolling is even greater, then T should be increased above 
.25 sec. Figure 8.11 compares the pitch response of the model 
to the data given by Cohen et. al. The above simulation was 
rerun with = 0 (representative of a real catapult launch) 

and the predicted perception of pitch was essentially the same. 


i 









Movement of DOWN for Catapult Launch Simulation 



Figure 8.11 Pitch Perception for Catapult Launch Simulation 
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8*3.4 Frequency Response for Small Pitch and Roll Oscillations 

The model 1 s use of otolith and canal information can 
be best understood by comparing the frequency response of the 
model to that of the sensors. Since the model has nonlinear- 
ities for large tilt angles and for conflicting sensory 
information, it is important to confine the oscillations to 
small angles (<10°) and to insure that only simple tilting 
or pitching stimuli are used. The response is essentially 
the same in both pitch and roll so only the data from the 
roll stimuli will be illustrated. Eight frequencies from 
.05 to 2.1 rad/sec were tested with stimulus amplitudes of 
5-10°. Lower frequencies were not tested since extremely 
long and therefore costly simulations would be necessary. 

Higher frequencies could not be tested since the update 
interval for the simulation was 1 second. Figure 8.12 shows 
the phase response. of the model for these frequencies. The 
amplitude response of the model is within 5% of unity over 
the range of frequencies tested. It is clear from Figure 
8.12 that for low frequency stimuli the model relies on oto- 
lith information and for higher frequency stimuli the model 
relies on information from the semicircular canals. The 
crossover frequency is approximately at .5 rad/sec. Nashner 
(Ref. 55 ) found a crossover frequency of approximately 
. 1HZ = .628 rad/sec from experiments involving postural 
control of pitch orientation. Since the phase and ampli- 




Figure 8- 12 Phase Response of Combined Model to Small Tilts (<icf) in Pitch and/or Roll 
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tude responses are so close to that of a unity gain for 
frequencies up to about 3 rad/sec the model predicts that 
our perception for small random tilt oscillations about a 
head erect position in Ig should be essentially correct. 
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8 • 4 Summary 

In this chapter a model for the perception of dynamic 
orientation resulting from stimuli which involve both the 
otoliths and the semicircular canals was developed. The model 
was applied both qualitatively and quantitatively to several 
such stimuli and its predictions evaluated. In all cases the 
model predictions were in substantial agreement with the 
known illusions or with the relevant experimental data. 
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CHAPTER XX 

CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 

Human perception of dynamic orientation based upon 
vestibular information has been modelled for near- threshold 
and suprathreshold stimuli. Before the processing of 
vestibular information could be modelled it was necessary to 
model the information available from the peripheral sensors. 
One conclusion which can be drawn from these models of 
afferent response is that while it is possible that little 
or no central processing is taking place for simple canal 
stimulation it is almost certain that significant dynamic 
processing is occurring in the case of stimuli which only 
involve otolith function. Furthermore , it has been demon- 
strated that the differences between- the afferent responses 
observed for the otoliths and the subjective responses 
seen in psychophysical experiments can be reconciled and that 
this reconciliation is consistent with the assumption of 
optimal processing by the higher centers. 

The conclusions which can be drawn from this research 
about the processing of vestibular afferent information by the 
brain are summarized in the following two sections. Finally, 
the chapter concludes with some suggestions for research 
which could extend the results presented in this thesis. 

9.1. Summary of Threshold Modelling 

Before a reasonable model could be constructed to predict 




the detection probabilities for arbitrary near threshold 
stimuli, it was necessary to determine the basic mechanism 
which gave rise to sensory thresholds. Two fundamentally 
different mechanisms were considered* The first hypothesis, 
called the "simple threshold model," consisted of a dead 
zone nonlinearity associated with the peripheral sensor which 
blocked the response from any stimulus which was not suf- 
ficiently large. The second hypothesis considered was that 
sensory thresholds arose only because the stimulus generated 
afferent response was masked by the variations in afferent 
firing which are independent of the stimulus* These hypotheses 
could be distinguished experimentally by determining the 
threshold level (75% correct detection) associated with a 
stimulus which is proportional to the sum of a subject’s 
velocity step and acceleration step thresholds. Such an 
experiment Was carried out and the results clearly demonstrated 
that the second hypothesis, designated the "signal in noise 
model," was to be preferred to the "simple threshold model.” 
Furthermore the data indicated that the phenomena of vestibular 
thresholds could be accounted for by a model of central 
processing of vestibular information consisting only of an 
optimal processing of afferent firing rates in additive noise 
with no necessity for peripheral dead zone nonlinear- 
ities . 

Once the signal in noise hypothesis was accepted as 
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an adequate model for the mechanism underlying the threshold 
phenomenon, it was necessary to develop a model of the 
processor which could predict the detection probabilities as 
a function of time for arbitrary near threshold stimuli. 

The model which resulted incorporates a first order 
processor which attempts to eliminate the tonic discharge 
and any associated low frequency, low amplitude variations 
so that the remaining signal can be assumed to consist of 
only the stimulus related signals and essentially uncorrelated 
measurement noise. This signal can then be processed 
sequentially to produce either a "moving right" response, a 
"moving left” response or no response at the end of each measure- 
ment interval. Through the use of Monte Carlo simulations 
with different sample functions of the noise process, a 
histogram of responses as a function of time can be generated 
for any given stimulus which should reveal both the total 
probability of correctly detecting that stimulus over a given 
period of time and the general distribution of response 
latencies. The threshold model for rotational stimuli 
correctly predicts the threshold magnitudes for velocity 
steps, acceleration steps and combination stimuli (velocity 
step plus acceleration step) and correctly indicates the general 

distribution of response latencies. Finally it is interesting 
to note that the detector used to model rotational stimuli 

when coupled with the afferent models for the otoliths was 
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found to be adequate for predicting both the acceleration 
step threshold for the utricles and the saccules and in 
addition gave reasonable predictions for the detection laten- 
cies at threshold . While it is possible that such a finding 
is coincidental it seems more likely that it is indicative 
of the fact that the higher centers detect motion by proces- 
sing near threshold canal and otolith information in a similar 
manner- 

9.2. Summary of Suprathreshold Modelling 

The problem of modeling human perception of supra- 
threshold stimuli was divided into three parts . The 
first part consisted of modelling the afferent information 
available from the sensors and coupling this with a model 
of central processing suitable for noninteracting stimuli. 

The results of this effort were threefold; 

1, Predictions could be made for the dynamic response 
to simple noninteracting stimuli , 

2- The best estimate of the head’ s rotational rate 

based upon information from the semicircular canals 
and the best estimate of the directionand magnitude 
of the specific force vector based upon otolith 
information was available for further integration 
for the case of interacting stimuli. 


and 
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3. A consistent mathematical framework had been 
developed for the central processor which incor- 
porated a model of the a priori information about 
the stimulus, a model of the sensory dynamics and a 
model of the variations in afferent firing, and which 
indicated that at least in the use of otolith informa- 
tion, the central processor made a significant 
contribution to the total dynamic response. 

The second part of this investigation centered on the 
perception of state orientation (no canal information) with 
respect to a constant specific force field. A thorough 
review of the illusions of static orientation indicated that 
they were consistent with a simple vector transformation which 
could be associated with differences in the processing of 
signals arising from stimuli in and stimuli perpendicular to 
the "utricle plane.” Based on these observations a model 
was developed which incorporated this difference in processing 
and which was capable of predicting the direction and magni- 
tude of the experimentally determined illusions of orientation. 
Finally it was observed that the alteration of saccular infor- 
mation required by the model was similar to the nonlinear 
response to static tilts seen by Fernandez et al. (Ref. 24) 
in otolith afferents in the squirrel monkey. This similarity 
suggests that the mechanism which gives rise to these illusions 
may have at least part of its origin in the peripheral sensor. 



Finally the problem of integrating i reformation from 
the semicircular canals and the otolith for the general 
class of interacting stimuli was considered* The major 
difficulty encountered in modelling the perception of dynamic 
orientation for motions which involve rotations about a hori- 
zontal axis was the problem of deriving the transformation 
of canal and otolith information which produces a perception 
of orientation with respect to the vertical* Once such a 
transformation is derived , predictions for the other per- 
petual outputs (rotation rates „ accelerations, etc.) follow in 
a relatively straightf orward manner* The model far the percep- 
tion of the vertical relies primarily on the otolith sensors 
for low frequency ( ^.5 rad/sec) changes in orientation 
and relies primarily on canal information which is confirmed 
by changes in the direction of the perceived specific force 
sensed by the otoliths for more rapid changes ( > "*~6 5 rad/sec) 

in orientation* This spectral division of responsibility 
is quite reasonable in light of the frequency characteristics 
of the sensors and the problems associated with any attempt 
to differentiate between translational accelerations and a 
change in orientation with respect to the vertical. The 
response of the model to small variations in tilt angle with 
respect to the vertical in a 1 g environment indicate that 
under these conditions the perception of tilt should be 
essentially correct for frequencies from zero to about three 
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radians per second. Accuracy in this region of operation 
should be expected since this is the region in which most 
head movements take place in daily life. 

The model's usefulness in predicting, without detailed 
simulation, the qualitative nature of the response to be 
expected from relatively simple interacting stimuli was 
demonstrated with several examples. Furthermore the accuracy 
of the model's quantitative predictions were shown for 
several stimuli for which data was available. While suggestions 
for further research to improve this model will be given in the 
next section, the results of the simulations carried out indi- 
cate that the model in its present form should be very useful 
in predicting the perceptual response to a wide variety of 
stimuli which up until now could not be confidently predicted. 

9.3. Suggestions for Further Research 

The results presented in this thesis could be extended 
by further research in the following areas: 

1. There is a great need for further information about 
the dynamic response of otolith afferents to time 
varying changes in specific force. The model of oto- 
lith information developed in Chapter Three is 
consistent with qualitative descriptions of otolith 
response but should be compared to more qualitative 
data. Specifically, a systematic study of the response 
of first order afferents to stimuli of various 
frequencies would be very useful. 
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2. Investigation of possible afferent vestibular 
thresholds might be conducted by studying the 
response of afferent fibers to stimuli which are 
near the perceptual threshold* One difficulty 
which would arise would be that the noise on a single 
afferent fiber would be significantly greater than 
the response due to the stimulus. The only way to 
circumvent this problem would be to average the 
responses of many stimulus trials. 

3. The statistics of afferent noise should be investi- 
gated much more thoroughly than has been done up to 
this point. Not only should the autocorrelation of 
the noise process on a single afferent be studied 
but also its correlation with the noise processes 
on other sensory afferents. Such an investigation 
should mdicate if the low frequency variation of the 
tonic discharge which was postulated in this thesis 
is present and would also suggest the degree to which 
afferent channels can be considered independent. 

4. Further psychophysical experiments should be carried 
out to determine the effect of a suprathreshold 
stimulus to one sensor on the detection probabilities 
associated with a near threshold stimulus to another 
sensor. The results of such a study might have 
significant implications concerning the applica- 
bility of subthreshold stimuli for the improvement of 
simulator fidelity. 
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5. The amount of quantitative data available from 
psychophysical experiments in which subjects are 
exposed to interacting stimuli and for which they 
are requested to indicate their perception of the 
vertical is quite limited. Any experimental program 
which systematically investigates these responses 
would provide very useful information for the model- 
ling of the perceptual response to interacting stimuli. 

6. Neurophysiologic studies of the interaction of 
semicircular canal and otolith information might 

be very productive. Since the otoliths provide the 
steady state response to continuous rotations about 
a horizontal axis in the same way that the visual 
system does for rotations about a vertical axis, 
it would be reasonable to assume that higher order 
vestibular neurons exist (most likely in the medial 
vestibular nucleus since this is the first nucleus 
which receives significant projections from both the 
semicircular canals and the otoliths) which depend 
upon both semicircular canal afferents and otolith 
afferents and which correctly reflect the true rotation 
about a horizontal axis {see Dischgans et al Ref. 22 arri 
Henn, Young, Finley Ref. 37 for evidence of such an 
interaction between semicircular canal and visual in- 
formation for rotations about a vertical axis) . 
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7 . Finally an attempt should be made to incorporate 

visual information into the model of human perception 
of dynamic orientation. While there is not sufficient 
information at this time to develop a definitive model 
of visual-vestibular perception of dynamic orientation 
there is a great need for a preliminary model which 
includes visual information 0 Such a model would be 
useful to suggest critical psychophysical and neurophysi- 
ological experiments and could be modified or if 
necessary radically altered in light of new experi- 
mental results. Even if the model underwent several 
radical changes the explicit nature of such a mathe- 
matical model helps to clarify both the issues in- 
volved and the underlying assumptions which too 
often are made but not explicitly recognized or under- 


stood. 
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Appendix A 


SUMMARY OF PARAMETERS FOR PERCEPTUAL MODELS 


I. SEMICIRCULAR CANAL MODEL 


The response of semicircular canal afferents is modelled 
by the following linear dynamical system (see section 3.1). 


*0 “ *<*5 + Sc" 


(A.l) 


and 


y c = 


5c 

03 


+ SFR + n 
c c 


where 


- 


§C - 


0 

0 

.37037 


-17.7966 -200.0888 


(A. 2) 
(A. 3) 

(A. 4) 

(A. 5) 
(A. 6) 


C = (-23.5785, -1131.89, -6371.86, 63.6620] 

SFR is the spontaneous firing rat© (90 ips) 
c including low frequency variations 

is a white measurement noise process which 
is discussed in the section 3.1 
u is the effective stimulus to the canal in (rad/sec) 
y c is the first order afferent response of the canal 


(ips) 

In tiie computer simulations used to test the model the 
afferent response was updated every .1 sec, by the following 
discrete model: 
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c (t+al) = ^ 


a) (t+.l) 


+ SFR + n 
c c 


(A . 8) 


where 


,9899912 
- . 17517x10 
— , 00183626 


09959244 .47296x10 


8915322 


00496791 


-.0879745 -.440444x10 


(A. 9) 


D c = 


2241507x10 

4730012x10' 

4958018x10' 


(A. 10) 


C , SFR and n are the same as previously defined and 
uj is the effective stimulus (rad/sec) 

Since the Kalman filter was designed to process one 
measurement each second and to estimate not only the internal 
states of the sensor but also the stimulus input, to , it 
must have a model of the canal dynamics and the input process 
This is accomplished by augmenting the state vector of equa- 
tion A.l with a state x^(4) which represents the input to (see 
figure 3.3). The resulting linear model which represents 
the internal model used by the Kalman estimator is given by: 


Sci “ ^ci^ci + 5 ci w ci 


(A. 11) 


y cx ^CI^CI * SFR CI + n CI 


(A. 12) 


-.37037 -17.7966 -200.0888 


(A. 13) 


0 


-5 
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B 


Cl 


(A. 14) 


and 


W CI is a unit white noise process (A- 15) 

y Cl' £ci' SFR CI / n CI are the same clS y c' Sc' 


SFR c' n c 


The discrete version of the Kalman filter (see equations 
3.11, 3.13, 3.14 and 3.15) requires the transition matrix 
associated with which is given by 


Cl 


.99911 . 

.95698 

.004760 

. 767xl0~ 3 

-.001763 

.91448 

.004573 

.924xl0~ 3 

-.001694 

-.07953 

- .. 415x10*" 3 

~.490xl0" 4 

0.0 

0.0 

o 

• 

o 

.006738 


(A. 16) 


The resultant steady state Kalman Gains are given by 

-3 


K. 


-.9191772x10 

+.8141562x10 

+.1503512x10 

+.3055998x10 


-5 

-3 

-1 


(A. 17) 


Finally the steady state Kalman filter for the semi- 
circular canal system is given by 


x(t+l) = + K fl9C [(y CI -SFE CI )-C <:i 0 CI x c (t)] 


(A. 18) 


where x^^t+l) is the minimum mean squared error estimate 

of the rotational rate u>(t+l) 
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y CI is the afferent measurement at t+1 


and 


^CI' 


-°°C f sF ^CI* ^CI are defined above 


IX 


OTOLITH MODEL 


The response of utricular afferents is modelled by 
the following linear dynamical system (see section 3.2) 


where 


and 


X,, 

= A x + 

B, f 


(A. 19) 

— u 

— u— u 




a C X + 
-u-u 

sm u 

+ n u 

(A. 20) 


\ ° 


1 1 



A » 
u 


B u K 


C u - 


-40. 


- 200.2 


[1300, 18000] 


(A. 21) 


(A. 22) 


(A. 23) 


SFR u is the spontaneous firing rate (88 ips) (A. 24) 


including low frequency variations 


n u is a white measurement noise process which 
is discussed in section 3.2 
f is the effective specific force acting on the 
sensor (g's) 

y u is the afferent response of the utricle (ips) . 
The discrete version of equations (A. 19) and (A. 20) 
for an update interval of .1 seconds is given by 


x u (t+.l) - $ u (t+.l,t)x u <t) + D u f (t+.li 


(A. 25) 


Y u (t+.l) = C u x(t+.l) + SFR u + n u 


(A. 26) 
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where 





. 98118 .0049059 

-.19624 -.00098118 

.47050xl0~ 3 
.49059x1 (f 2 


(A. 27) 
(A. 28) 


C u , SFR u and n u are as previously defined and f is 
the specific force stimulus (g's). The internal model of 
the otolith dynamics augmented with the internal model of 
the stimulus statistics is given by 


*01 ^JI-UI + 

y UI := ^UI^UI + SFR UI + n UX 


where 



0 1 0 

-40. -200. 2 1 

0 0 -1 



(A. 29) 
{A. 30) 

(A. 31) 
(A. 32) 


W yi is a unit. white noise process (A. 33) 

and y UJ , C UI , SFR UIf n u;[ are the same as y y , C 
SFR^j and hy 


The transition matrix associated with A u;[ for a 1 


second update interval is given by: 
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.81955 . 0040977 

- 0 16391 -.81955x10 

0.0 QoO 


. 0028114 
* 0012864 
. 35788 


(A. 34) 


The steady state Kalman gains for the system defined by 
(A. 29) and (A. 30) are 


K 


U 


.1283x10 

.2792x10 


-4 

-2 


(A. 35) 


.6102x10 

Finally the steady state Kalman filter for the utricle 
is given by 


— U ( t + ^ ) £| + ^ ^UI -SFR UI^ ~ "UI— UI^U (A. 36) 


y yi is the afferent measurement at t+1 
and $ 0I , K^, Sfr ui# , are defined above. 

The saccular model is identical to that of the utricle 

except that the afferent response is only half as great (y = 

s 

y^/2) and the resulting Kalman gains are twice as great (K^- 
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